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An  analytical  model  was  developed  to  describe  the  thermodynamic 
and  fluid  dynamic  processes  in  an  exhaust  gas  cooler  employing  liquid 
water  injection.  The  model  is  based  on  the  solution  of  the  equations  of 
conservation  of  species,  momentum,  and  energy  for  the  system  and  the 
equations  for  the  exchange  of  these  quantities  between  liquid  and 
gaseous  phases.  These  equations  are  programmed  for  solution  on  an 
IBM  360  computer.  The  predictions  of  the  model  are  compared  with  mea¬ 
sured  data  from  a  series  of  turbojet  tests  in  the  Propulsion  Development 
Test  Cell  (T-l)  spray  cooler.  The  comparison  showed  that  the  model  gave 
a  good  agreement  with  the  measured  pressure  and  liquid  temperature  at 
various  points  along  the  cooler.  Parameters  such  as  gas  temperature  and 
specific  humidity  which  were  not  measured  are  discussed  in  terms  of 
their  relation  to  the  overall  cooler  performance.  From  the  results  of 
the  measurements  and  predictions,  a  physical  description  of  the  cooling 
process  is  presented.  Based  on  the  results  of  one  of  the  tests,  a 
possible  method  of  reducing  the  pressure  loss  in  a  cooler  is  proposed. 
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ABSTRACT 


An  analytical  model  was  developed  to  describe  the  thermodynamic 
and  fluid  dynamic  processes  in  an  exhaust  gas  cooler  employing  liquid 
water  injection.  The  model  is  based  on  the  solution  of  the  equations  of 
conservation  of  species,  momentum,  and  energy  for  the  system  and 
the  equations  for  the  exchange  of  these  quantities  between  liquid  and 
gaseous  phases.  These  equations  are  programmed  for  solution  on  an 
IBM  360  computer.  The  predictions  of  the  model  are  compared  with 
measured  data  from  a  series  of  turbojet  tests  in  the  Propulsion  Develop¬ 
ment  Test  Cell  (T-l)  spray  cooler.  The  comparison  showed  that  the 
model  gave  a  good  agreement  with  the  measured  pressure  and  liquid 
temperature  at  various  points  along  the  cooler.  Parameters  such  as 
gas  temperature  and  specific  humidity  which  were  not  measured  are 
discussed  in  terms  of  their  relation  to  the  overall  cooler  performance. 
From  the  results  of  the  measurements  and  predictions,  a  physical  de¬ 
scription  of  the  cooling  process  is  presented.  Based  on  the  results  of 
one  of  the  tests,  a  possible  method  of  reducing  the  pressure  loss  in  a 
cooler  is  proposed. 
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SECTION  I 
INTRODUCTION 


Testing  of  turbojet  engines  and  rocket  motors  at  simulated  altitude 
in  ground  test  facilities  requires  cooling  of  the  high  temperature  ex¬ 
haust  gas  to  a  relatively  low  temperature  before  the  gas  enters  the  ex¬ 
haust  gas  pumping  system.  Cooling  of  the  gases  by  water  spray  with 
direct  heat  and  mass  exchange  between  the  water  and  the  exhaust  gas 
has  been  utilized  in  many  test  facilities.  This  method  of  cooling  is 
often  called  spray  cooling. 

Many  of  the  spray  coolers  used  in  the  Engine  Test  Facility  (ETF) 
at  the  Arnold  Engineering  Development  Center  (AEDC)  receive  exhaust 
gas  from  a  rocket  or  turbojet  engine.  The  cooling  process  reduces  the 
temperature  from  approximately  4000°R  (maximum  temperature  of  a 
turbojet  engine  exhaust  gas)  to  approximately  550°R.  By  means  of  an 
atomizing  water  spray,  the  exhaust  gas  is  cooled  and  humidified.  The 
cooling  produces  a  temperature  compatible  with  the  ducting,  control 
valves,  and  pumping  system  material  limits.  Water  conservation  is 
an  important  consideration  in  operation  because  of  the  large  quantities 
of  spray  water  required. 

Previous  work  has  developed  computer  models  for  spray  coolers 
based  on  the  assumption  of  a  homogeneous  two-phase  flow  with  kinetic 
and  thermodynamic  equilibrium  (Refs.  1,  2,  and  3).  The  investigations 
contained  in  this  report  cover  the  development  of  a  computer  model  of 
a  spray  cooling  process  that  follows  a  typical  liquid  water  droplet  in  the 
cocker  ducting.  No  assumptions  of  kinetic  or  thermodynamic  equilib¬ 
rium  between  the  gas  and  liquid  are  made.  The  model  is  then  compared 
with  measurements  made  in  a  spray  cooler  during  operation.  The  approach 
approach  is  similar  to  that  used  by  Shapiro  (Ref.  4). 


SECTION  II 

EXHAUST  GAS  COOLING  SYSTEM 


2.1  CONFIGURATION 

The  configuration  of  the  Propulsion  Development  Test  Cell  (T-l) 
spray  cooler  consists  of  a  diverging  conical  inlet  section  followed  by 
a  constant-area  duct  to  the  end  of  the  spray  cooler  (Fig.  la. 

Appendix  I).  The  cooling  water  is  introduced  through  a  group  of  nozzles 
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arranged  in  a  series  of  banks,  in  which  the  first  three  banks  of  sprays 
consist  of  nozzles  projecting  a  fan-type  spray  directed  downstream 
along  the  wall  to  protect  the  ducting  (Fig.  lb).  The  remaining  banks 
are  arranged  in  a  wagon-wheel  configuration  with  several  spokes,  each 
"spoke"  containing  several  spray  heads.  Each  spray  head  contains 
several  fixed-geometry,  conical  spray  nozzles  (Figs,  lb  and  c)  directed 
generally  downstream.  The  water  to  each  spray  bank  is  supplied  by  a 
large  header,  and  the  flow  rate  to  each  spray  bank  is  controlled  by  a 
valve  between  the  header  and  the  spray  bank. 


2.2  INSTRUMENTATION 

Instrumentation  was  provided  to  measure  flow  rates  and  pressures 
of  the  exhaust  gas  stream  entering  the  cooler,  and  e  temperature  and 
composition  were  calculated  using  the  method  of  R.  ..  5.  Exhaust  gas 
static  pressure  and  liquid  water  temperature  measurements  were  auso 
made  at  five  axial  stations  along  the  cooler.  The  temperature  of  the 
cooling  water  before  injection  was  measured,  and  the  flow  rate  was 
calculated  from  pressure  measurements  made  across  an  orifice  or  at 
the  control  valve.  The  location  of  this  instrumentation  is  shown  in 
Fig.  2,  Measurements  taken  by  this  instrumentation  provided  experi¬ 
mental  correlation  with  the  analytical  results  from  the  mathematical 
model. 

The  millivolt  outputs  from  the  thermocouples  and  strain-gage-type 
pressure  transducers  were  recorded  on  either  magnetic  tape  or  by  a 
photographically  recording  galvanometer-type  oscillograph.  The  mag¬ 
netic  tape  data  were  reduced  on  a  digital  computer,  and  the  oscillograph 
data  were  reduced  manually  using  electrical  calibrations  taken  prior  to 
testing. 
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SECTION  III 

DEVELOPMENT  OF  THE  ANALYTICAL  MODEL 


The  analytical  model  is  baseu  .  (  the  equations  of  conservation  of 
energy,  momentum,  and  species  for  the  ekhaust  gas  cooler  and  the  ex¬ 
change  of  these  quantities  between  phases.  The  model  considers  the 
behavior  of  a  typical  drop  down  the  length  of  the  copier  and  calculates 
the  changes  in  thermodynamic  properties  over  a  series  of  small  incre¬ 
mental  distances.  The  equations  were  programmed  for  solution  on  a 
digitial  computer. 

Some  of  the  key  assumptions  in  the  analysis  are  as  follows: 

1.  All  gases  including  water  vapor  obey  the  perfect  gas 
equation  of  state. 

2.  The  flow  is  steady  and  one  dimensional. 

3.  The  gat.  mixture  at  any  section  is  homogeneous. 

4.  The  droplets  from  each  injection  station  are  uniformly 
distributed  over  the  cross-sectional  area  of  the  cooler. 

5.  The  droplets  are  injected  parallel  to  the  ghs  flow  and 
maintain  this  direction  throughout  the  cooler.  The  in¬ 
fluence  of  gravity  on  the  droplets  is  considered  negli¬ 
gible. 

6.  There  is  no  aerodynamic  breakup  or  agglomeration  of 
the  drops. 

7.  The  drops  injected  at  any  injection  station  are  uniform 
in  size. 

8.  The  maximum  number  of  injection  stations  is  nine, 
and  their  spacing  is  arbitrary. 

9.  The  internal  resistance  of  the  drops  to  heat  distribu¬ 
tion  is  negligible,  thus  the  temperature  is  uniform 
through  the  drop. 

10.  The  drops  injected  at  each  injection  station  are 
accounted  for  separately. 

11.  Heat  transfer  and  friction  at  the  duct  walls  and  piping 
are  negligible. 

12.  Cross-seeliunal  area  of  the  cooler  is  a  prescribed  func¬ 
tion  of  distance  along  the  cooler. 
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3.1  EQUATION  FOR  THE  CONSERVATION  OF  SPECIES  FOR  ONE 
INJECTION  STATION 

The  conservation  of  species  written  for  the  exhaust  gas  and  water 
(in  both  liquid  and  vapor  form)  over  an  incremental  distance  fix  is 

dm  dm.  dm 

m  +  rh.  +  m  =  m  +  — 3— -  fix  +  m.  +  -3 — fix  +  m  +  — -r -c-  fix  ^ 
v  £  11c  v  dx  i  dx  nc  dx 

1  • 

If  the  noncondensable  flow  rate  is  assumed  constant  and  insoluble  in 
water,  Eq.  (1)  may  be  simplified: 

dm  dm„ 

_£+_J=0  (2) 

dx  dx 

The  mass  fraction  of  vapor  may  be  written 


C..  = 


v  m  +  m 

v  nc 


and  the  mass  fraction  of  noncondensable  gas  is 


1  q  = _ _ _  =  1  _  c 

nc  m  +  m  v 

v  nc  . 

In  addition,  the  specific  humidity  may  be  defined  as 


1  -  C 


and  the  liquid  water  ratio  (f^)  as 


"  m 


3.2  EQUATION  FOR  THE  CONSERVATION  OF  ENERGY  FOR  ONE 
INJECTION  STATION 

The  energy  equation  is  developed  to  equate  the  total  energy  of  the 
exhaust  gas  and  liquid  watdr  as  they  pass  two  planes  an  incremental 
distance  (fix)  apart. 
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m  [h  +(V  2/2)l  +  m  fh  +  (V  2/2)l  +  m,  [h  +  (V.2/2)l 
ncLnc  nc  J  v  L  v  nc  J  i  (_  *  *  J 

=  [m  +  (dm  /dx)6xlih  +  (dh  /dx)6x 
[  nc  nc  J  l  nc  nc 

+  [Vnc  +  (dVnc/dx)6x]2  /2}  +  K  +  (drfiv/dx)6x] 

+  <dh  /dx)6x  +  Tv  +  (dV  /dx)6xl2/2l 
v  L  nc  nc  J  ; 

+  |rn^  +  {dm^/dx)6x|  jh^  +  (dh^/dx)6x 
+  jV£  +  (dV£/dx)6x]2/2| 


By  expanding  and  simplifying  the  above  and  expressing  the  mass  flow 
rates  of  the  various  components  in  terms  of  Eqs.  (5)  and  (6),  Eq.  (7) 
becomes 


i  dV  C  dh  dV 


x  nc  dx  1  -  C  dx  nc  dx 

v 


h  +  (V  12)  dC 


(1  -  c/ 


(^♦v  ♦  !*-)**.  . 

\dx  £  dx  J  X  £  2  /  dx 


3.3  EQUATION  FOR  THE  CONSERVATION  OF  MOMENTUM  FOR  ONE 
INJECTION  STATION 

The  momentum  equation  expressing  the  total  momentum  passing 
two  planes  (6x)  apart  is: 

m  V  +  m  V  +  m.  V  +  PA  =  [rh  +  (dm  /dx)6xl  fv 
nc  nc  v  nc  £  £  L  nc  nc  J  L  nc 

+  (dV  /dx)6xl  +  l"mv  +  (dm  /dx)6x]fv  +  (dV  /dx)6x 
nc  J  L  v  v  J  L  nc  nc 

+  |rh^  +  (dm^/dx)6xj  +  (dV^/dx)6xj 
+  [p  +  (dP/dx)6xJ  ^A  t  (dA/dx)6x] 
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Expanding  and  simplifying  give 

rft  (dV  /dx)  +  V  (dm  /dx)  +  m  (dV  / dx)  +  V  (dm  /dx)  (10) 
nc  nc  nc  nc  v  nc  nc  v 

+  mi(dVjf/dx)  +  V^drf^/dx)  +  P(dA/dx)  +  A(dP/dx)  «  0 


By  dividing  by  and  incorporating  the  specific  humidity  and  liquid 

water  ratio  in  terms  of  Eqs.  (5)  and  (6),  respectively,  Eq.  (10)  may  be 
expressed  as 


dV  C  dV 

nc  ,  v  nc 


nc 


dC 


dx 


1  -  C  dx 
v 


(i  -  c  y 

v 


dV,  df„  _ 

v  .  ,  _ i  +  v  — i.  +  P  dA 

dx  *  dx  i  dx  m  dx 

nc 


,  A  dP 

m  dx 
nc 


(ID 


Multiplying  by  V  (1  -  C_)'  and  assuming  that  the  change  in  area  over 
the  increment  6x  is  negligible  give 


dV 


dC 


(1  -  C  )V 


,  .  +  V  :  — j —  +  (1  -  C  )'  V 

v  nc  dx  nc  dx  v  nc 


If  Hi  +  v  —1 

dx  s.  dxj 


+  a'Cy>  >  o 

p  dx 

nc 


(12) 


3.4  EQUATIONS  FOR  MULTIPLE  INJECTION  STATIONS 


Exhaust  gas  coolers  similar  to  those  shown  in  Fig.  la  consist 
of  a  series  of  spray  banks  or  water  injection  stations,  whereas  the 
equations  previously  developed  indicate  that  the  water  is  injected  uni¬ 
formly  at  one  station.  The  equations  arc  easily  expanded  to  include 
multiple  injection  stations  by  adding  a  term  to  describe  the  liquid  injec¬ 
tion  conditions  at  each  spray  bank  and  the  location  of  each  bank.  The 
previously  developed  equations  (Eqs.  (2),  (8),  and  (12))  are  modified 
as  shown  below.  Equation  (2)  becomes 


dm 


dm{ 

dx 


=  0 


(13) 
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Equation  (8)  becomes 


and  Eq.  (12)  becomes 


(1  -  C  )  V 

v  nc  dx 


dV  dC 

+  V  — -  +  (1  -  C  )  V  £  f 
dx  nc  dx  v  nc  “ 


- 1  +  V 

i .  ox  £. 

l 


dti] 

.  dx 
i 


(1  -  C  )  ^ 

- —  -o 

p  dx 


The  piping  necessary  for  the  spray  banks  in  an  exhaust  gas  cooler 
similar  to  the  one  shown  in  Fig.  la  can  occupy  a  significant  portion  of 
the  cross-sectional  area  of  the  duct.  In  the  cooler  of  test  cell  T-l,  the 
frontal  area  of  the  piping  at  each  injection  station  is  approximately  12 
percent  of  the  total  cross-sectional  area.  Because  of  the  large  amount 
of  liquid  normally  used  for  cooling  and  the  blockage  caused  by  water 
piping,  the  equations  describing  the  cooling  process  must  be  further 
modified  to  incorporate  terms  necessary  to  account  for  the  loss  in 
momentum  of  the  liquid  that  strikes  this  piping.  This  modification  is 
made  0.25  ft  upstream  of  each  spray  bank  where  the  liquid  properties 
from  all  previous  stations  are  mass  averaged  to  produce  two  new 
streams,  one  of  which  represents  the  liquid  passing  a  spray  bank  with¬ 
out  interference  and  a  second  stream  which  strikes  the  piping,  losses 
its  momentum,  and  is  then  reaccelerated.  The  fraction  of  liquid  strik¬ 
ing  the  piping  is  equal  to  the  fraction  of  area  occupied  by  the  piping. 
Therefore,  the  new  liquid  properties  passing  the  spray  bank  will  be 


4  Y,  Y +  'Y 


(10) 

(17) 
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Ts  V 

t 

n  1 


(18) 


D|  f.  +  D>  f.  +  Dj  f .  + 


da  = 


L  lt. 

n  i 


(IP) 


The  properties  of  the  liquid  that  strikes  the  piping  will  be 

V  -1.0 

B 


T  =  T 
SB  SA 


13a 


p  V  -  V, 
g\  nc 


(20) 

(21) 

(22) 


The  velocity  (V  )  is  arbitrarily  set  at  a  small  value  but  not  zero  be- 
B 

cause  it  appears  in  the  denominator  of  several  calculations;  T„  is 

B 

assumed  equal  to  T  since  the  system  is  adiabatic  and  no  heat  is  lost 
bA 

to  the  piping.  The  diameter  of  the  reaccelerated  drop  (Dg)  is  based  on 

Eq.  (12.  6)  of  Ref.  6.  The  equation  has  been  modified  by  neglecting  the 
second  term  on  the  right  side  since  it  is  almost  negligible  for  the  condi¬ 
tions  encountered  in  this  program.  The  final  form  of  the  equation  for 
the  drop  diameter  is  basically  a  solution  to  the  Weber  number  for  a 
critical  value  of  13. 


3.5  EQUATIONS  FOR  THE  EXCHANGE  OF  MASS,  ENERGY,  AND  MOMENTUM 
BETWEEN  PHASES 

It  is  now  necessary  to  develop  the  equations  to  relate  the  transfer 
of  mass,  energy,  and  momentum  between  phases.  Since  the  mass  flow 
rate  of  the  noncondensable  portion  of  the  exhaust  stream  is  considered 
constant  and  negligibly  soluble  in  water,  the  only  exchange  of  masc 
occurs  between  the  liquid  water  and  vapor.  The  conservation  of  water 
was  expressed  earlier  in  Eq.  (2)  and  is 
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dm 
_ v 

dx 


+  E 


dm 


dx 


=  0 


Furthermore  the  mass  transfer  to  or  from  a  single  drop  may  be 
developed  from  Eq.  (21.2-26)  of  Ref.  7  which  expresses  the  molar  rate 
of  evaporation  as 


=  k  ttD.2 
x.  1 
i 


1  -  X 

vs. 

1 


(23) 


By  multiplying  by  the  molecular  weight  of  vapor  (M  ),  the  results  in 
terms  of  the  mass  rate  of  evaporation  may  be  expressed  as 


u 


(m) 


v. 

l 


M 

v 


X  -  X 
vs.  v 
l 


jtD.*  M  ■ 
l  v  1  -  x 


vs . 


(24) 


Since  the  mass  rate  of  evaporation  is  equal  to  the  decrease  in  the  mass 
of  a  drop  per  unit  time,  then 


dM , 
d 

dt 


-  k  nD.  M 
x.  1  v 
i 


(25) 


Since  the  equations  developed  will  be  solved  for  incremental  distances 
(dx),  the  equation  above  will  be  more  useful  in  terms  of  the  distance  (dx) 


=  k  ttD.2 
x.  1 


M 


x  -  x 
V  vs. 

1 


v  1  -  x 


vs. 

1 


(26) 


For  a  mixture  of  perfect  gases,  the  molar  concentration  (x  )  can  be 

i 

expressed  in  terms  of  the  pressure  where 

P 

vs. 

X  =  -5-i  (27) 

vs.  P 
1 
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is  the  mole  fraction  of  vapor  at  the  drop  surface  and  Pvg  is  the  vapor 

saturation  pressure  computed  at  the  drop  surface  temperature.  While 
this  method  of  evaluating  xvg  is  exact  only  for  zero  mass  transfer,  it 

can  be  shown  to  give  satisfactory  results  even  at  relatively  high  mass 
transfer  rates.  The  mole  fraction  (xv>  of  the  free  stream  is 


x 

V 


m  /M 
v  v 


C  /M 

V  V 


(m  /M  )  +  (m /M  j 
v  v  nc  nc 


'M  )  +  (1  -  C  /M  ) 
v  v  nc 


(28) 


The  change  in  the  total  amount  of  liquid  may  be  expressed  in  terms  of 
the  change  for  one  drop  and  the  number  of  drops 


"V  dMd. 

_ l  _  _ i_  _ i 

dx  M ,  dx 
a. 
l 


(29) 


The  transfer  of  energy  between  phases  is  related  to  the  thermody¬ 
namic  state  of  the  exhaust  gas  stream  and  the  liquid  drops.  Since  the 
system  is  adiabatic  and  at  a  constant  area  over  the  distance  (dx),  any 
change  in  the  gas  stream  will  necessarily  result  in  a  change  in  the 
drops;  therefore,  the  exchange  of  energy  between  phases  will  be 
expressed  as  a  change  in  internal  energy  of  a  single  drop,  and  this  will 
then  be  related  to  the  change  in  energy  of  all  the  liquid.  The  change  in 
internal  energy  for  a  single  drop  over  the  distance  (dx)  is 


where  the  first  term  on  the  right  expresses  the  convective  heat  transfer 
and  the  second  term  expresses  the  heat  transfer  accompanying  the  mass 
transfer  and  change  in  phase.  From  the  known  energy  transfer  for  one 
drop,  the  total  energy  transfer  may  be  expressed  as 
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It  is  assumed  that  the  resistance  to  heat  distribution  within  the  droplet 
is  negligible  compared  with  the  resistance  to  heat  transfer  at  the  sur¬ 
face,  that  is,  the  temperature  within  the  drop  is  uniform.  Thus, 


c*(T  -  T  J 
i  s.  ref 


(32) 


for  a  constant  specific  heat  of  liquid  (c^). 


The  momentum  transfer  between  phases  will  be  expressed  using 
Newton's  Second  Law  where  the  force  on  the  drop  is  due  only  to  the 
droplet  drag.  Therefore,  expressed  in  this  manner, 


which  may  be  simplified  to 


(34) 


where  p  is  the  density  of  the  combined  noncondensable  gas  and  vapor 
& 

in  the  stream. 


The  equation  of  state  for  the  exhaust  gas  stream  is 


where 


(p  +  p  )R  T 
nc  v  g  nc 


(35) 


(36) 


3.6  COMPUTER  SOLUTION  OF  THE  EQUATIONS 

The  conditions  in  the  cooler  are  determined  by  first  computing  the 
changes  to  the  liquid  and  then  incorporating  these  changes  into  the  solu¬ 
tion  of  the  conservation  equations  for  the  complete  system.  The 
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computer  solution  is  based  on  the  modified  Euler  method.  The  changes 
in  the  liquid  properties  are  calculated  from  the  derivatives  given  in 
Eqs.  (26),  (29),  (30),  and  (34)  and  a  known  step  size  (dx).  These 
changes  are  calculated  for  each  liquid  stream.  The  sum  of  the  changes 
in  mass,  energy,  and  momentum  are  then  incorporated  into  Eqs.  (13), 
(14),  (15),  and  (35)  to  solve  for  the  new  gas  properties.  An  iteration 
technique  is  used  for  the  solution  to  the  last  three  eouations.  The  calcu¬ 
lation  procedure  continues  down  the  cooler  until  a  point  is  reached 
0.  25  ft  upstream  of  a  spray  bank.  At  this  point,  the  liquid  properties 
are  averaged  as  discussed  in  Section  3.4  (Eqs.  (16)  through  (22)).  Bv 
using  the  gas  properties  last  calculated  and  the  liquid  properties  of 
stream  "a"  only,  the  calculation  procedure  discussed  above  is  com¬ 
pleted  for  one  step  (dx).  At  this  point,  stream  "b"  (the  liquid  that  has 
impinged  on  the  piping)  is  added  to  the  calculation,  and  the  changes  to 
Eqs.  (26),  (29),  (30),  and  (34)  are  calculated  for  the  two  streams 
separately.  These  changes  in  mass  energy  and  momentum  are  included 
in  this  iterative  solution  to  Eqs.  (13),  (14),  (15),  and  (35).  The  calcula¬ 
tion  procedures  continue  until  a  new  spray  bank  is  reached  and  then 
these  liquid  properties  are  included  in  the  calculation  routine. 

The  step  size  (dx)  is  variable  in  this  program.  The  initial  value 
used  is  0.0001  ft,  but  if  convergence  is  achieved  quickly  (less  than  3 
iterations),  the  step  size  is  increased  for  the  next  series  of  calculations. 
The  step  size  will  vary  between  0.  0001  and  0.  01  ft  depending  on  the 
number  of  iterations  necessary  for  convergence  in  the  previous  set  of 
calculations.  A  computer  listing  of  the  program  is  given  in  Appendix  III. 

The  frequency  of  printout  for  the  calculations  may  also  be  varied, 
but  experience  has  shown  that  for  most  conditions  printing  the  results 
every  0.25  ft  is  sufficient  to  see  the  changes  in  the  exhaust  gas  cooler 
conditions. 

Typical  input  for  a  computer  run  is  shown  in  Tables  la  and  b 
(Appendix  II).  Special  note  should  be  taken  of  the  spray  banks  in  which 
no  water  is  injected  (f^  =  0).  The  spray  banks  are  included  in  the  input 
because  they  will  contribute  blockage  to  the  system  and  their  location 
must  be  known.  In  each  case,  a  fictitious  velocity,  temperature,  and 
drop  size  is  also  included  to  prevent  division  by  zero  during  the  solu¬ 
tion,  but  since  these  properties  are  also  multiplied  by  f  ,  thev  become 
zero  and  do  not  affect  the  final  solution. 

For  conditions  where  two-phase  flow  exists,  but  the  piping  for  in¬ 
jecting  the  liquid  does  not  occupy  a  significant  portion  of  the  cross- 
sectional  area,  a  variation  of  the  analytical  model  may  be  used.  This 
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variation  involves  only  the  solution  of  Eqs.  (13),  (14),  (15),  (2G),  (29), 
(30),  (34),  and  (35)  without  the  averaging  of  the  liquid  properties  dis¬ 
cussed  in  Section  3.  4  (using  Eqs.  (16)  through  (22)).  In  addition,  a 
drop  size  distribution  may  be  simulated  in  this  variation  of  the  model 
by  inputing  the  various  drop  diameters  and  their  respective  quantities 
(f .  )  as  spray  stations  but  with  the  stations  at  the  maximum  dx  distance 
i 

apart.  A  computer  listing  for  this  model  variation  will  be  found  in 
Appendix  IV. 


SECTION  IV 

EVALUATION  OF  THE  ANALYTICAL  MODEL 


An  evaluation  of  the  computer  model  was  made  by  comparing  the 
predicted  exhaust  gas  pressure  with  measured  data  and  also  the  pre¬ 
dicted  liquid  water  temperature  with  the  value  measured  by  an  exposed 
junction  thermocouple  at  several  points  in  the  exhaust  gas  cooler.  In 
addition,  the  calculated  exhaust  gas  temperature  and  the  specific  humid¬ 
ity  are  discussed  to  determine  how  these  parameters  are  influenced  by 
test  conditions. 

Six  typical  data  points  taken  during  the  testing  of  a  turbojet  engine 
are  used  to  evaluate  the  model.  The  input  conditions  for  use  in  the  com¬ 
puter  program  are  shown  in  Tables  I  and  11.  Also  shown  in  the  tables 
are  the  inlet  conditions  to  tin*  spray  cooler  which  were  calculated  using 
the  method  of  Ref.  4  which  lias  been  included  as  a  part  of  the  model. 
These  data  show  that,  for  the  runs  in  Table  I,  the  cooler  inlet  conditions 
are  constant  and  the  difference  is  in  the  spray  banks  that  are  in  use, 
whereas  the  data  in  Table  II  show  the  cooling  water  parameters  to  be 
constant  and  the  exhaust  gas  temperature  and  velocity  to  vary.  Two 
other  items  of  importance  that  should  be  noted  are: 

1.  The  cooling  water  flow  rate  from  the  individual  "wagon- 
wheel  ’  spray  banks  was  kept  constant,  and  only  the  num¬ 
ber  of  spray  banks  was  varied,  and 

2.  The  cooling  water  from  the  wall  sprays  of  spray  banks 
No.  2  and  3  were  not  normally  included  in  the  calculation, 
whereas  the  water  from  spray  bank  No.  1  was  included. 

The  flow  rate  from  the  individual  spray  banks  was  kept  constant  to  mini¬ 
mize  i  lie  effects  of  variations  iii  cooling  w  ater  velocity,  drop  size,  and 
distribution  from  the  individual  nozzles.  The  cooling  water  from  spray 
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banks  2  and  3  was  not  included  because  it  was  believed  that  these  wall 
sprays  would  not  contribute  significantly  to  the  cooling  of  the  exhaust 
gas.  Calculations  show  that  the  pressure  change  in  the  ducting  to  the 
first  wagon-wheel  can  generally  be  adequately  described  by  assuming 
a  one -dimensional  isentropic  flow  with  no  cooling  water  present.  The 
cooling  water  from  spray  bank  No.  1  was  included  because  liquid  water 
must  be  present  at  the  start  of  the  computer  program  (i.  e. ,  f^  0)» 

and  the  difference  between  the  isentropic  value  of  pressure  and  that 
calculated  using  the  initial  spray  bank  was  not  significant. 

The  measured  and  calculated  values  of  pressure  as  a  function  of 
distance  along  the  cooler  are  shown  in  Fig.  3;  the  liquid  water  tempera¬ 
tures  are  shown  in  Fig.  4.  The  increase  in  pressure  during  the  first 
9  ft  (the  diverging  portion  of  the  ducting)  followed  by  a  drop  in  pressure 
for  the  constant  area  portion  of  the  cooler  is  characteristic  of  nearly 
all  runs.  The  initial  rise  in  pressure  is  due  primarily  to  the  subsonic 
compression  of  the  exhaust  gas  with  very  little  acceleration  of  the  cool¬ 
ing  water  except  on  the  outer  edges  of  flow.  The  abrupt  decrease  in 
pressure  that  follows  occurs  in  the  constant  diameter  section  of  the 
cooler  where  the  wagon-wheel  sprays  are  located.  The  drop  in  pres¬ 
sure  in  this  section  is  caused  by  acceleration  of  the  cooling  water  from 
the  initial  wagon-wheel  spray  bank  and  the  loss  in  momentum  of  the 
previously  injected  liquid  water  as  it  strikes  the  piping  and  is  then  re¬ 
accelerated.  This  later  loss  in  momentum  is  taken  in  account  by  the 
averaging  of  the  liquid  properties  and  the  resultant  use  of  Eqs.  (34), 

(35),  and  (36).  Although  the  area  occupied  by  the  internal  water  piping 
at  each  spray  station  is  small  (approximately  12  percent),  it  is  sufficient 
to  cause  a  drastic  change  in  the  pressure  characteristics.  The  magnitude 
of  the  change  caused  by  inclusion  of  this  piping  is  best  illustrated  by 
assuming  that  in  the  constant  diameter  section  the  piping  is  removed  but 
the  water  is  still  introduced  uniformily  over  the  cross  section  of  the 
cooler  at  the  various  spray  banks.  Figure  5  shows  the  calculated 
cooler  pressure  for  several  blockages  as  well  as  the  standard  12  percent 
used  for  the  calculation  of  all  data  from  the  cooler  in  test  cell  T-l.  All 
calculated  pressures  show  good  agreement  initially,  but  then  the  pres¬ 
sure  begins  to  increase  for  the  case  of  zero  blockage  while  decreasing 
rapidly  for  the  remaining  cases.  This  increase  in  static  pressure  is 
caused  by  the  decrease  in  dynamic  pressure  while  the  total  pressure  of 
the  exhaust  stream  remains  essentially  constant.  The  decrease  in 
dynamic  pressure  is  due  primarily  to  the  cooling  of  the  exhaust  gas 
stream.  The  decrease  in  pressure  for  the  conditions  with  various 
amounts  of  blockage  is  due  to  the*  interference  caused  by  the  piping. 

The  effect  of  blockage  in  a  flow  stream  is  well  known  and  the  above 
example  illustrates  its  importance  in  a  two-phase  stream.  The  percent¬ 
ages  in  Fig.  5  cover  the  range  of  normal  and  extreme  blockage  conditions 
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and  indicate  not  only  the  importance  of  including  the  blockage  in  the  ana- 
lytical  model  but  also  the  importance  of  minimizing  it  wherever  possible 
in  cooler  design. 

The  measured  and  predicted  liquid  temperatures  for  the  data  points 
in  Table  I  are  shown  in  Fig.  4.  Three  of  the  four  measured  values 
show  good  agreement  with  the  predicted  values,  while  the  remaining 
value  (the  initial  measurement)  is  always  high.  It  is  significant  to  note 
that  the  predicted  liquid  temperature  for  the  first  wagon-wheel  spray 
bank  at  9.  2  ft  rises  approximately  70  °R  in  approximately  6  in.  and  then 
levels  off  at  an  almost  constant  value  even  when  additional  spray  banks 
are  used.  This  rapid  rise  in  temperature  is  due  to  the  fact  that  initially 
the  liquid  water  temperature  is  low  (536°R)  and  its  vapor  pressure  at 
the  drop  surface  is  also  low.  The  low  vapor  pressure  of  the  drop  com¬ 
bined  with  the  low  vapor  partial  pressure  in  the  gas  stream  results  in 
a  low  mass  transfer  rate,  while  the  large  difference  in  gas  and  liquid 
temperatures  gives  a  high  heat  transfer  rate  and  a  rapid  rise  in  the 
temperature  of  the  liquid  with  little  evaporation.  As  the  liquid  temper¬ 
ature  begins  to  rise,  the  difference  between  the  partial  pressure  of  the 
drop  and  gas  stream  increases,  and  the  mass  transfer  (or  evaporation 
from  the  drop)  increases.  This  increase  in  mass  transfer  continues 
until  a  liquid  temperature  is  approached  where  heat  transfer  to  the  drop 
is  almost  completely  used  for  the  evaporation  of  water.  The  final  tem¬ 
perature  approached  by  the  liquid  is  its  adiabatic  or  wet  bulb  saturation 
temperature.  As  additional  cooling  water  is  added  through  the  use  of 
additional  spray  banks  this  process  is  repeated,  but  the  rate  of  liquid 
temperature  rise  will  decrease  because  of  the  smaller  temperature 
difference  between  gas  and  liquid  and  also  because  the  hotter  liquid 
must  also  be  cooled.  The  cooling  of  the  liquid  does  not  become  im¬ 
portant  until  there  is  a  very  large  amount  of  "hot"  liquid  present.  The 
fact  that  the  cooler  liquid  temperature  does  rise  very  rapidly  keeps  the 
overall  cooling  process  from  becoming  extremely  inefficient  due  to 
alternately  heating  and  cooling  the  liquid  in  the  stream. 

The  previously  mentioned  thermocouple  located  at  9.  2  ft  always 
reads  high.  The  high  reading  is  believed  to  be  caused  by  the  location 
of  the  thermocouple  at  the  edge  of  the  diverging  section  where  the  fan- 
type  spray  will  leave  a  liquid  deficient  region  near  the  thermocouple. 
Since  the  smaller  liquid  quantity  is  surrounded  by  a  large  amount  of 
hot  exhaust  gas,  the  heat  transfer  to  the  liquid  will  be  abnormally  high 
and  thus  cause  the  liquid  temperature  to  rise  to  a  value  higher  than  is 
predicted  by  the  computer  model  which  assumes  a  uniform  liquid  distri¬ 
bution  over  the  cross  section  of  the  cooler. 
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The  validity  of  the  computer  model  is  best  determined  by  comparing 
the  predicted  and  measured  values  of  static  pressure,  gas  temperature, 
and  specific  humidity.  A  comparison  using  the  static  pressure  has 
already  been  made,  but  measurements  of  the  later  two  quantities  have 
not  been  made  because  of  the  difficulties  inherent  in  a  two-phase  stream. 
The  exhaust  gas  temperature  is  important  because  of  the  effect  on  pump¬ 
ing  machinery  capabilities.  As  the  temperature  of  the  exhaust  gas 
entering  the  machines  is  increased,  the  maximum  mass  flow  rate  that 
can  be  pumped  at  a  constant  pressure  decreases;  or,  stated  another  way, 
for  a  given  mass  flow  rate  of  exhaust  gas,  the  minimum  upstream  pres¬ 
sure  increases  as  the  temperature  of  the  exhaust  gas  entering  the 
machinery  increases.  Therefore,  it  is  desirable  to  cool  the  exhaust 
gas  as  much  as  possible.  The  specific  humidity,  like  the  gas  tempera¬ 
ture,  places  a  lower  limit  on  the  pressure  capabilities  of  the  exhaust 
machinery.  As  the  specific  humidity  increases,  the  minimum  upstream 
pressure  at  the  test  cell  also  increases  because  this  additional  vapor  is 
additional  mass  that  must  be  removed.  Therefore,  the  optimum  condi¬ 
tion  would  appear  to  consist  of  the  lowest  exhaust  gas  temperature  and 
specific  humidity.  The  problem  is  that  the  temperature  normally  de¬ 
creases  at  the  expense  of  an  increase  in  humidity  for  spray  coolers  like 
those  in  FTF  unless  very  large  quantities  of  water  are  injected.  Since 
the  overall  process  of  reducing  the  temperature  is  normally  by  evapora¬ 
tion  of  cooling  water,  the  specific  humidity  for  the  process  increases  as 
the  temperature  decreases. 

The  predicted  gas  temperatures  along  the  length  of  the  cooler  is 
shown  in  Fig.  6  for  the  three  data  points  under  discussion.  The  de¬ 
crease  in  temperature  follows  the  same  path  for  the  three  runs  as  long 
as  the  same  spray  banks  are  used.  As  expected  the  lowest  gas  temper¬ 
ature  occurs  for  the  data  point  using  the  most  cooling  water  (Run  No. 
36-13),  while  the  highest  temperature  is  predicted  with  the  least  amount 
of  cooling  water  (Run  No.  36-13).  The  temperature  curve  shows  a  dis¬ 
tinct  change  occurring  at  approximately  9  ft.  Although  the  data  indicated 
that  the  pressure  change  at  this  point  could  be  treated  as  a  subsonic  com¬ 
pression  of  a  gas  with  no  mass  transfer,  this  is  not  the  case  with  the 
predicted  cooling  curve.  If  the  gas  temperature  followed  a  subsonic 
compression  process  with  no  mass  transfer,  the  predicted  temperature 
should  be  at  some  value  greater  Ilian  the  cooler  inlet  value  of  approxi¬ 
mately  3570°R.  If  this  were  a  subsonic  compression  the  wall-type 
sprays  would  probably  be  spraying  directly  along  the  wall  and  the 
thermocouple  at  9.  2  ft  should  be  reading  approximately  gas  temperature 
and  the  first  wagon-wheel  spray  bank  should  be  surrounded  by  the  high 
temperature  gas  flow.  As  noted  earlier,  the  thermocouple  at  the  end 
of  the  divergent  section  indicates  a  measured  temperature  higher  than 
the  predicted  liquid  but  certainly  not  an  exhaust  gas  temperature.  Since 
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the  thermocouple  is  measuring  a  liquid  temperature  (Refs.  8  and  10), 
there  has  obviously  been  heating  of  the  water  indicating  that  the  temper¬ 
ature  characteristics  cannot  be  described  by  an  isentropic  compression. 
Therefore,  some  cooling  and  mass  transfer  has  taken  place  in  the  diver¬ 
gent  section  of  the  cooler,  and  the  predicted  curve  probably  has  the 
correct  shape,  but  it  is  not  possible  to  know  if  the  temperature  is  abso¬ 
lutely  correct.  The  remainder  of  the  predicted  exhaust  gas  tempera¬ 
ture  curve  is  typical  — a  rapid  decrease  in  temperature  while  there  is  a 
large  temperature  difference  between  gas  and  liquid  followed  by  a  der 
creasing  rate  of  cooling  near  the  exit  as  the  temperature  difference  de¬ 
creases.  The  point  where  the  cooling  curves  separate  is  the  location 
of  the  next  spray  bank  being  used.  '  The  differences  noted  4t  the  exit 
indicate  the  magnitude  of  temperature  change  that  can  be  expected  by 
using  additional  spray  banks  for  the  conditions  of  these  tests. 

Since  one  of  the  objects  of  the  cooling  process  is  to  get  the  lowest 
value  of  exhaust  gas  temperature  at  the  lowest  specific  humidity,  the 
temperature  as  a  function  of  the  specific  humidity  is  presented  for  the 
three  data  points  in  Fig.  7  to  show  how  the  cooling  takes  place.  The 
two  points  immediately  obvious  and  important  to  the  model  description 
are: 

1.  The  overall  process  is  one  of  humidification  and  not  two 
separate  processes,  i.  e.  ,  one  of  humidification  followed 
by  dehumidification,  and 

2.  There  are  short  periods  of  dehumidification  downstream 
of  each  active  injection  station  but  quickly  followed  by  a 
resumption  of  the  evaporative  process. 

The  normal  method  of  visualizing  the  cooling  process  in  an  exhaust 
gas  spray  cooler  is  to  picture  first  a  short  section  of  cooler  in  which 
sufficient  water  is  present  to  provide  saturation  conditions.  This  water 
is  injected  into  the  stream  and  is  immediately  vaporized  because  of  the 
large  temperature  difference  between  liquid  and  exhaust  gas.  This  pro¬ 
cess  involves  transferring  sufficient  heat  from  the  gas  stream  to  vapor¬ 
ize  the  water.  After  the  gas  stream  becomes  saturated  with  respect  to 
the  cooling  water  that  has  been  injected,  any  further  cooling  water 
serves  to  dehumidify  the  gas  stream.  This  dehumidification  process  is 
generally  imagined  to  take  place  very  slowly  when  compared  with  the 
vaporization  process.  As  shown  in 'Fig.  7,  the  cooling  process  does  not 
appear  to  follow  the  model  described  above;  instead  the  process  appears 
to  be  one  of  almost  continuous  evaporation  with  a  few  periods  of  slight 
dehumidification.  Which  process  is  correct  becomes  very  important 
in  the  understanding  and  design  of  spray  coolers.  The  first  process 
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actually  describes  what  takes  place  in  an  infinitely  long  cooler  where  1 
6nly  small  amounts  of  iwater  are  injected,  and  this  water  is  allowed  to 
reach  temperature  and  velocity  equilibrium  before  any  more  water  is 
injected.  ’  In  this  way,  saturation  m4y  be  achieved  but  With  no  excess  ,, 
liquid  water  present.  The  second  process  describes  a  nonequilibritim 
process  in  terms  of  liquid  and  vapor  temperatures^  velocity,  and  con¬ 
centrations  but  is  the  abtual  process  in  an  exhaust  gas.  cooler. 

,  When  liquid  water  is  injected  into  the  gas  stream,  initially  the 
liquid,  is  at  a  low  temperature  and  partial  pressure,  whereas  the  gas  i 
stream  has  a  relatively  high  temperature  but  low  vapor  pressure  due 
to  t^ie  small' amount  of  vapor  in  the  stream  (generally  only  the  water 
formed  during  the1  combustion  process).  $ince  the  evaporation  pro-  . 
cess  is  controlled  by  the  vapor  concentration  difference  ,(see  Eq.  (16)),  ' 
the  mpss  transfer  will  initially  be  very  low,,  but  because  of  the  large 
temperature  differences  (on  the  order  of  30d0°R),'  the  heat  transfer 
rate  will  be  very  high.  With  the  high  heat  transfer  rate,  the  liquid 
temperature  r,ises  very  rapidly,  and  the  partial  pressure  difference 
between  the  liquid  and  gad  stream  increases,  causing  a  rise  in  the 
mass  transfer  rate.  Thie  process  of  rising  liquid  temperature  and 
'mass  transfer  rate  continues  until  a  liquid  temperature  is  approached 
where  the  heat  transferred  intoithe  liquid  is  used  almost  corripletely  for,  ■ 
vaporization.  The  temperature  that  is  approached  by  the  liquid  is  the 
adiabatic  or  wiet-bulb  saturation  temperature,  but  although  this  temper¬ 
ature  is  nearly  achieved,  the  evaporation  of  the  liquid  continues  because 
there  'still  exists  a  partial  pressure  difference  between  the  droplet  sur¬ 
face  and  gqs  stream  tp  provide  the  mass  transfer  driving  force;  and  a 
temperature  difference  to  supply  heat  for  vaporization.  This  process 
will ’continue  until  the  partial  pressure1  and  temperature  difference  dis¬ 
appears. 

This  process  dis'eussed  above  is  basically  for^ne  spray  bank,  where 
as  the  normal  cooler  operation  usels  several  banks.  What  happens  when 
fresh  cooling  water  (from  a  downstream  spray  bank)  is  injected  into  the 
gas  stream  can  be  divided  into  two  processes.  These  occur:  , 

'  '  j  *  t 

1.  When  the  vapor  pressure  of  the  fresh  liquid  is  greater 
than  the  vapor  pressure  of  the  exhaust  stream,  and 

2.  Wljen  the  vapor  pressure  of  the  fresh  liquid  is  less  , 
than  the^  vapor  pressure  of  the  gas  streapi. 

I  ; 

i  Both  of  these  conditions  are  shown  in  Fig.  7c.  The  first  condition 
occurs  at  the  entrance  to  the  cool  dr ’when  the  exhaust  gas  contains  very 
little  vapor  (the  partial  pressure  is  almost  negligible)  and  water  is 


1  ° 

-L  u 


AEOC-TR72-89 


injected  with  a  partial  pressure  of  approximately  64  psf.  In  this  instance, 
evaporation  begins  immediately  and  continues  the  length  of  the  cooler. 

The  second  case  occurs  when  the  gas  temperature  has  reached  2500aR 
(at  spray  bank  No.  4).  The  partial  pressure  of  the  gas  stream  is  now 
215  psf  which  is  above  that  of  the  incoming  water  and  should  result  in 
the  gas  stream  being  dehumidified.  This  is  in  fact  what  happens,  but 
because  the  dehumidification  takes  place  over  such  a  short  distance,  the 
decrease  in  specific  humidity  does  not  show  up.  The  dehumidification 
process  is  much  clearer  at  spray  bank  No.  5  when  the  gas  temperature 
has  reached  1600°R,  and  the  specific  humidity  decreases  from  0.  43  to 
0.41  before  beginning  to  increase  again.  Another  way  of  picturing  the 
process  is  shown  in  Fig.  8  where  the  partial  pressure  difference  be¬ 
tween  the  liquid  and  vapor  is  shown  as  a  function  of  the  specific  humid¬ 
ity.  When  the  pressure  difference  is  positive,  evaporation  takes  place, 
and  when  it  is  negative,  dehumidification  takes  place.  In  this  figure, 
the  liquid  from  spray  bank  No.  1  is  shown  to  be  evaporating  from  the 
start,  whereas  the  liquid  from  spray  bank  No.  4  initially  causes 
dehumidification  until  the  pressure  difference  reaches  zero.  Then  the 
process  for  spray  bank  No.  4  becomes  evaporative,  and  the  specific 
humidity  begins  to  increase  again.  The  same  thing  will  happen  to  the 
spray  banks  downstream  as  shown  in  Fig.  8  where  dehumidification 
takes  place  immediately  downstream  of  the  injection  station.  The  spe¬ 
cific  humidity  has  a  greater  decrease  for  each  succeeding  spray  bank 
because  the  partial  pressure  difference  is  initially  greater. 

Three  additional  runs  are  included  in  this  discussion  to  show  the 
agreement  between  the  model  and  the  experimental  data  and  also  to 
point  out  some  possible  effects  of  the  two  wall  spray  banks  in  the  diverg¬ 
ing  section  which  are  normally  omitted  from  the  calculation  procedure. 
The  measured  engine  inlet  parameters  and  the  calculated  cooler  inlet 
parameters  are  shown  in  Table  II,  and  the  pressure  as  a  function  of 
distance  data  is  shown  in  Fig.  9.  These  data  show  good  agreement 
between  theory  and  experiment  at  the  cooler  exit,  but  for  Run  3G-1G,  the 
agreement  at  the  exit  from  the  diverging  section  (9.  2  ft)  is  not.  good. 
Comparing  the  measured  and  predicted  pressure  for  the  three  runs 
shows  that  the  agreement  is  good  for  Run  3G-14,  but  becomes  progres¬ 
sively  poorer  for  Runs  3G-15  and  3G-1G.  From  Table  II,  the  test 
parameters,  including  the  cooling  water,  are  seen  to  be  the  same  with 
only  the  engine  fuel  flow  rate  changing.  Thus,  there  is  poorer  agree¬ 
ment  between  the  model  and  data  as  the  fuel  flow  rate  decreases. 

A  possible  explanation  for  this  lies  in  the  interaction  between  wall 
spray  banks  2  and  3  and  the  relatively  low  \-c1ocity  gas  stream.  At  the 
interface  between  the  liquid  and  gas  streams,  a  portion  of  the  hot  ex¬ 
haust  gas  is  cooled  by  flowing  radially  through  the  wall  sprays  to  the 
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area  betweer  the  interface  and  the  duct  wall  while  the  remainder  under¬ 
goes  essentially  an  isentropic  compression  as  it  flows  down  the  duct. 

As  the  Mach  numbers  of  the  two  streams  decrease,  their  static  pres¬ 
sure  increases,  and  for  the  gas  flowing  through  the  sprays,  the  total 
pressure  will  also  increase  because  of  the  evaporation  of  water.  The 
resultant  effect  is  to  increase  the  static  pressure  above  that  predicted 
by  the  model  because  of  the  cooling  of  the  radially  flowing  gas. 

To  show  the  effect  of  the  wall  spray  banks  on  the  agreement  between 
measured  and  predicted  pressure,  a  run  (35-44,  Table  II)  was  chosen 
with  similar  cooler  inlet  conditions  but  without  spray  banks  2  and  3 
operating.  The  predicted  and  measured  pressure  for  the  first  9.  2  ft  is 
shown  in  Fig.  10  for  Runs  35-44  and  3G-1G.  The  agreement  between 
predicted  and  measured  pressure  for  Run  35-44  is  good,  indicating  that 
the  wall  sprays  are  at  least  part  of  the  problem. 

A  second  interesting  point  about  Run  3G-16  is  that  the  measured 
pressure  is  higher  than  the  total  pressure  calculated  by  assuming 
either  an  isentropic  expansion  or  a  two-phase  cooling  process,  as 
shown  in  Fig.  11.  To  achieve  this  measured  pressure,  mass  must  be 
added  to  the  stream  but  under  conditions  where  there  is  not  significant 
loss  in  momentum  due  to  the  added  mass  being  accelerated.  This  con¬ 
dition  could  probably  be  achieved  by  the  gases  flowing  radially  through 
the  water  stream  and  adding  vapor  to  the  gas  stream. 

The  condition  of  abnormally  high  static  pressure  is  not  usually 
encountered  because  the  wall  sprays  arc  not  used  for  conditions  such 
as  encountered  in  Rim  3G-1G  where  the  cooler  inlet  velocity  and  temper¬ 
ature  are  low.  For  those  runs  wiiere  the  sprays  are  used,  such  as 
3G-14,  the  total  energy  level  of  the  stream  masks  any  influence  of  the 
wall  sprays. 

If  this  radial  mass  flow  through  the  wall  sprays  is  the  reason  for 
the  extremely  high  static  pressure  at  the  end  of  the  diverging  section, 
this  technique  might  be  a  useful  way  of  supplementing  the  exhaust  pump¬ 
ing  machinery  to  achieve  a  lower  test  cell  pressure. 


SECTION  V 

CONCLUDING  REMARKS 


An  analytical  model  was  developed  to  describe  the  process 
on  by  an  exhaust  gas  spray  cooler.  The  model  consisted  of  the 
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solution  to  tlic'  equations  of  conservation  of  species,  momentum,  and 
energy  and  the  exchange  of  these  quantities  between  the  gas  and  liquid 
phase. 

The  model  was  compared  with  data  from  turbojet  tests  conducted  in 
Propulsion  Development  Test  Cell  (T**!).  The  range  of  spray  cooler 
inlet  conditions  was  as  follows: 

m  =  148  to  155  lbm/sec 
nc 

T  =  10G6  to  ■J5ou°R 
nc 

P,  =  711  to  009  psf 

V  =  402  to  1147  ft/sec 
nc 

f  =  2.  5  to  3.  3 
£ 

T  =  536°R 
s 

The  static  pressure  and  liquid  temperature  were  measured  at  six  posi¬ 
tions  along  the  length  of  the  cooler,  and  these  pressures  and  tempera¬ 
tures  were  cempared  with  similar  values  predicted  by  the  model.  The 
measured  static  pressure  at  the  cooler  exit  agreed  witli  the  predicted 
value  within  4  percent  or  less,  and  the  measured  liquid  temperature 
agreed  within  1  percent  of  the  predicted  value.  With  the  exception  of 
the  pressure  at  the  entrance  to  the  cylindrical  section,  the  agreement 
between  model  and  measurement  for  the  other  data  was  at  least  this 
good.  The  predicted  values  of  gas  temperature  and  specific  humidity 
are  discussed,  but  measured  values  of  these  quantities  wei’e  not  avail¬ 
able  for  comparison  because  of  the  lack  of  adequate  instrumentation. 
Successful  instrumentation  was  not  available  for  making  this  type  of 
measurements  in  a  typical  exhaust  gas  cooler  stream  with  liquid-to-gas 
mass  ratios  on  the  order  of  2  to  1  or  greater. 

Static  pressure  measurements  at  the  exit  to  the  diverging  section 
of  the  cooler  gave  abnormally  high  results  for  one  run  based  on  the  pre¬ 
dicted  value  from  the  model.  A  possible  explanation  for  these  data 
based  on  a  separated  recirculating  flow  in  this  area  was  postulated. 
Additional  data  for  verifying  this  were  not  available.  A  possible  method 
for  improving  spray  cooler  performance  based  on  this  postulate  was 
also  mentioned. 
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Fig.  3  Predicted  and  Measured  Exhaust  Gas  Cooler  Pressure  as  a  Function 
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a.  Run  No.  36*13 

Fig.  4  Predicted  and  Measured  Liquid  Temperature  as  a  Function  of  Axial  Location 
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c.  Run  No.  36-18 
Fig.  4  Concluded 
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8.  Run  No.  36-13 

Fig.  7  Predicted  Relation  between  Exhaust  Temperature  and  Specific  Humidity 
for  Three  Cooling  Water  Flow  Rates 
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Fig.  8  Effect  of  the  Difference  in  the  Partial  Pressures  of  thfe  Liquid  and 
Exhaust  Gas  Streams  on  the  Specific  Humidity  for  Run  No.  36-13 
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Comparison  of  the  Measured  Cooler  Pressure  with  Several  Values  of  Calculated  Pressure 


TABLE  I  (Concluded) 
b.  Cooling  Water  Flow  Rate  Conditions 


Inlet  Liquid 


Run 

No. 

Station 

or 

Spray  Bank 

Distance 

from 

Entrance, 
x.  ft 

Liquid 
Ratio,  ff, 
lbm/lbm 

Velocity 
Vi.  ’ 
ft/ sec 

36-13 

1 

0 

0. 

524 

■n 

4 

9.25 

0. 

849 

mm 

5 

11.50 

890 

MM 

G 

13.  75 

0.  951 

95.0 

7 

1G.00 

0.0 

90.0 

C 

18.25 

0.  968 

94.0 

9 

20.  50 

0.0 

90.0 

10 

22.  75 

0. 

0 

90 

.0 

11 

25.00 

0. 

3 

90 

.0 

36-14 

1 

0 

0.520 

4 

9.25 

0.849 

91 

.  5 

5 

11.  50 

0.  908 

91 

.0 

6 

13.  75 

0.  9G1 

95 

.  5 

7 

16.  00 

8 

18.  25 

9 

20.  50 

10 

22.  75 

11 

25.00 

36-18 

1 

0 

0.  53 1 

71 

.0 

4 

9.  25 

0.  887 

92 

.0 

5 

11.  50 

0.  935 

91 

.  5 

6 

13.  75 

0 

7 

16.00 

8 

18.25 

9 

20.  50 

10 

22.  75 

. 

11 

25.00 

Temperature 


1  s> 


53  G 


Dri 


0.  21 
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36-15  4.32  147.23  I  779.0  5.58 


I 


I 


1 


TABLE  II  (Concluded)  * 
b.  Cooling  Water  flow  date  Conditions 


AEDC-TR-72-89 


Inlet  Liquid 


Station 


Spray  Bank 


Distance 
frofm! 
Entrance, 
x,  ft' 


0.0 
4.8 
,  G.  9 

9.25  I 

11.5 

,  13.75 
16..  00 
18.25 
20.  5 

0.0 
4'.  8 
G.  9 
9.  25 
11.  5 
13.  75 
1G.  0 

18.25 

20.5 


Liquid 
Ratio,  ff. 
Ibm  /  lWi 


0.  52 
0.  187 
0. 176  : 
0.849 
0.  908 
0,  961 
0.0 
0.0 
0.0 

0:534 
0.  191 
0.  179 
0.873 
0;  921 
0.  979 
0.0 
0.0 
0.0 

0.547 
0.  195 
0.  18G 
0.  910 
0.  959 
1.01 
0.0 
0.0 
0.0 

0.  175 
0.0 
0.  0 
1.  020 
1.341 
n.n 


I 


Velocity  Temperature,  Drop  Size, 


Vi. 

ft/ sec 


D.  ft' 


0.0020 
0.0020  , 
0.0020 
0.00059 


0.00197 
0.00197 
0.  00197 
0.  00059 


0.00194 
0.00194 
0:  00194 
0.00059 


0.0020 
0.  0020 
0.0020 
0.00059 
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APPENDIX  III 

A  LISTING  OF  THE  COMPUTER  PROGRAM 
AND  THE  REQUIRED  AUXILIARY  EQUATIONS 
FOR  AN  EXHAUST  GAS  COOLER 


A  listing  of  the  computer  program  for  the  solution  of  the  equations 
developed  in  the  text  is  given  in  this  section.  In  addition,  auxiliary 
equations  necessary  to  define  certain  constants  used  for  the  computer 
solution  are  listed. 


AUXILIARY  EQUATIONS 


These  auxiliary  equations  are  used  to  define  certain  dimensionless 
numbers  which  are  in  turn  solved  for  certain  coefficients  used  in  the 
equations  developed  in  the  text.  The  properties  of  the  system  used  in  the 
equations  are  evaluated  at  the  so-called  "film  temperature"  which  is  de¬ 
fined  as 

T  +  T 

T  _ S. 

fi  2 


where  the  f  indicates  a  film  property  and  the  i  identifies  the  particular 
liquid  stream  being  discussed.  For  this  program,  the  noncondensable 
gas  was  assumed  to  be  dry  air,  and  the  various  constants  were  calcu¬ 
lated  using  the  properties  of  air.  The  desired  constants  and  equations 
are  given  below: 
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C  /M 

-  V  V 

XV  =  C  /M  +  (1-  C  )M 

v  v  v  nc 


x 

V  . 

SI 


p 

V  . 

sx 

p 


V 


p 


fi 


(p)(Mfi) 
R  T„. 


COMPUTER  LISTING 

The  following  computer  program  was  programmed  for  the  IBM  360 
computer  and  was  used  to  obtain  the  calculated  data  of  this  program. 


I 


i 


\ 

1 

\ 
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DETERMINATION  OF  EXHAUST  CAS  COOLER 
INLET  CONDITIONS 

IMPLICIT  REAL*8<A-H,0-Z) 

REAL**  T2THI150) 

REAL**  AREA,RAIR,RFUEL 
REAL**  artm,arpr,arwi 
REAL**  I  RUN ( 2  >  * WHAT  II*) 

- CUNWOW/mOT  AR'IwnSUTtAMPtn  rrUTTTnTWTTTSUTTWPri - 

COMMON  t  ENTH  /  T.  TT,  R,  C0A(7>,  H(6),  C02AI7),  H2AI7),  XN2AI7I, 

1  02 A I  7 ) ,  H20 A ( 7 ) ,  WT<6>,  C0BI7),  C02B<7>,  H2B(7),  XN2SI7), 

2  02M7I,  H20B<7>,  CI6),  A(6),  P<212> 

CALL  ERRSETI261, 256,-1,1) 

WT<1  >  *  28.01 1 

- gyyyj  -  4gt01l - 

MTI3)  •  2.016 

MTI4)  ■  28.016 
WT<5)  ■  32.0 
WTI6)  s  IB. 016 
PI201 )■  11.766 

- P(2(ra>-,;  177?5T5 - 

PI205)*  12.770 
P( 207  >■  13.298 
P < 209  > ■  13.8** 

P(2ll  )*  1*  .*08 
Re  1.98726 

Z - PftPS5i)ftC  RATA — — TT5 — HTP - 

C 

RE AO  <5,100)  4P< I >,1*32,199) 

RE AO  <5,100)  < P< I ) , I >200,212 ,2  > 

100  FORMAT  18010.0) 

C 

Z  fEMPERATURE  COePF  ftl ENTS 

C 

READ  15.101)  ICOAl  l),Iel,7),  (C0B1 I >t 1*1,7),  ( C02 A ( I ) , I • 1 «  7  ) , 

1  (C02B1I ), 1  =  1.7),  I H2 A < 17,1*1.7),  < H2B ( I ),I»1»7)»  <  XN2 A4 I >,1*1,7), 

2  <XN2B<|),|«1,7»,  <02 A< I ) , I • 1 , 7 ) ,  < 02B1I) , 1 » 1 . 7 ) ,  <  H20A ( I ) , I • l , 7 ) , 

3  <H20R<I), 1-1,7) 

161  FORMAT  (5016.7/2016.7) 

C 

C  INPUT  CONDITIONS 

987  FORMAT  (20A*) 

C 

C  MASS  FRACTIONS 

Z - 

105  FORMAT  (6012.5) 

l  RE AO  <17,gN0«999>  C ,P 1 , T 1 , TW 1 , Vfil , VW 1 , TGUESS , I RUN, WHAT 
READ  (17)  AREA.RAIR.RFUEL 
TGUESSeTW l*.0l*(Tl-Twi ) 

T2MAX«TM1-.01 

- r2mN«TWIv.7M - - - 

IF  <P1 .EO.O.ANO.Tl.EO.O. )  GO  TO  999 

- WPtTffi - 

C  _ _ 

5M*0. 

00  6  I  «  l  ,6 

EM  =  EM  +  C  <  I  >/  MT( I > 
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CMC  1  •  1.0  -  CVl 
PVt  *(Fl  *  EM  *  C V 1  > /  1 A . 
PNC  1  «■  PI  -  PVl 
— fwkic - « i p i  ■>  pm  *  ccn  / 
T*T 1 /  l.R 
CALL  ENTHAL 
SUM  »  0.0 
DO  T  I  *  1.  5 
SUM  »  SUM  ♦  c (  I  >  *  H(  I  ) 

7~ COWTTNUc - 

HNCl  =  SUM  /  CNCl 


O0TM*RAIR«-RFUEL 

CONSTN=( ( 1545.*00TM)/( 28.85*144. *P1*AREA))**2 
CONSTN*CONSTN/(2.*77R.*32.2) 

TGS1*T1/1 .A 


'  u  jc  -  \  i  i-v/«  i  i  i  r  i  «n 

T  =  TGS 1 
CALL  ENTHAL 

CALL  SUMIT  (C.H.CNCl.ENTl ) 
T*T* 1  .A 

FUNC1SENT1  *CONSTN*T  *T-HNc'  1 


CALL  ENTHAL 

CALL  SUMIT  (C.H.CNCI.ENT2) 
T=T*1.R 

FUNC  2*6  N  T  2 ♦CONST  N*  T*  T  — MNC l 
401  TG|F*(TGSl*TGS2)/2. 


TsTGIF 

- CALL  ENTHAL - 

CALL  SUMJT (C.H.CNC1.ENT3) 

T  =  T*  1  •  R 

FUNC3*ENT3-,C0NSTN*T*T-HNC1 
TSTA  =  FUNC  l  *F(|mC3 
T$TR=FUNC2*FUNC3 

- tT"  itstai  404,4U4,402 - - 

402  IF  ( TST8 )  405,405,403 

403  WRITE  (6,411) 

411  FORMAT! *0  NO  ROOT  FOR  INITIAL  ENThALPM') 
GO  TO  1 

404  FUNC2  =  FUNC  3 

- rns2sTr,tF - 

GO  TO  401 

405  FUNC 1 =FUNC 3 
TGS1=TGIF 
GO  TO  401 

410  ENTNU=ENT3 

- rm^TTF  *ttp - 

C  END  OF  ITERATION 

VGI  =  (D0TM*1545.*T1 ) / ( P 1  * AR E A  *  1 44.  ) 

VG1 =VG1 /28 . 85 

WRITE  (6,413)  ENTNU.TI.VG1 

WRITE (9)C (6) ,VG1 ,T1 ,P1 

413  POftMAT  I  •  0 — EnThAiW  =  •  ,E12.4,  • - ?T5T_T 

*  •  GAi  VELOCITY  *  • ,E12.4) 

GO  TO  l 

999  CONT INUE 
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E NOFILE  0 
SfFPTWT” 
RETURN 
- 6N'D“ 


- gUHlOUTWETirRrC - 

IMPLICIT  REAl*8IA-H,0-Z > 

COMMON  7  ENTH  f  TT"TT  »  R  ♦  COA|7>,  HI6),  C02AI7),  H2AI7),  XN2AJ7), 

1  02AI7),  H20A ( 7  >  *  WT 1 6 )  ,  C08I7),  C02SI7),  H2BI7),  XN2BI7), 

2  02RI7  )  r  "M2081 7  f «  "CT6  V,  A(6>,PI212> 

3  TT -R*T *1*8 

- Dr-rr^jrorCT-iT)  to - 

00  4  J«I*6 

~ TTAT  J  J-C1TATT5 - 

CALL  hTRT  ( T ,  A,H( | ) ) 

'  HII)«HU>*TT/WTIU  *2848873 — *T7g7Frm - ~ - 

00  $  J-1,6 

5  A(~J>  '■  C02AI JV  ' 

CALL  HTRT  <T,A,HI2>> 

- HT?TaHT?  riTT/HT ( g ;  ♦yg?<sgv  *r.ff/WTlTI - 

00  6  J»l,6 

6  Af JI-H2AIJ) 

CALL  HTRT  (T.A.H13) ) 

- "HT3 ) «Hl3)*TT/wT(3)*2023,8-*l . 87WTI3 1 

00  7  J-1,6 

— r ‘anr-nreAm - 

CALL  HTRT  IT  , A , Hi  4  ) ) 

H(4)»Hm*TT/WT(4>  +2072,3  *1.8/WTI41 — — . 

00  8  J-1,6 

B  AIJ1-02AIJT  . . . .  . . . 

CAU  HTRT  IT,A,H<5>> 

- Hp>)-^(^TT|1TtfT|b,  -t?074.  - 

00  9  J-1,6 
9  AlJ)=H20Arj> 

CALL  HTRT  I T , A , HI  6  > ) 

H 1 6 ) -HIS1  *TT7WTT5T~*60164';T  "*1. 8  AWTC6) 

HV2-HI6) 

- GTT-nrrT - 

10  no  li  j-i,6 

11  A I J )  -COB  ( JT  .  . .  . 

CALL  HTRT  IT, A, Hlln 

HI  1  > -HI  1  >*TT/HT { l )  *28488.3"  M.8/WTU) 

00  12  J-1,6 

12  Al  JJ*<,u2«I  J)  ' 

CALL  HTRT  IT,A,HI2>) 

H|2)-HI?)*TT/WTI2)  *96290.  *I.8/WTTTr - 

00  13  J-1,6 

13  Al  JJ-H2BI J>  - 

CALL  HTRT  I T , A ,HI 3 ) ) 

- HI3  J  in 1 3 )  *TT /ITT  HT  1-202 i". S' "^lTR/W  1(3) - 

00  14  j-1,6 

14  Al JJ-XN2BI J) 

CALL  HTRT  (T,A,HI4>> 

H|4)=H|4)*TT/HTI4)  +2072,3  •1.B/WTI4) 

00  15  J-1,6 

15  M3l?n?R|JI - - - - - 

CALL  HTRT  ( T, 4, HIS) I 

HI5)=H(5)«tT/WTI5)  *2074.7  M.8/WTI5) 
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00  16  J»l,6 

15  AIJ)*h206(J)  — — . 

CAU  HTrT  (  T  ,  A , h( 6  )  ) 

rr r / w r t a i  huha.t  *i.»/wmbi 

17  CONTINUE 
HNC2  *0.0 
00  1*  t»l,5 

IB  HNC2»HNC2  ♦  C(I)*H(|) 
WNC,<«HNt2/TT.U~L16n 
RETURN 

- IHTO - 


SUBROUTINE  SUM IT (C*H,CN»ENTh> 

TNPi.  lull  REAL *B  (  - 

0 1  HENS  I  ON  C ( 1 > (HI  1 > 

S(IM*0. 

00  1  1*1(5 

1  SUN»SUM*Ct I)*HI I) 

ENTMaSUM/CN 

- RETURN - - -  — - 

END 
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MAIN  PROGRAM  FOR  EXHAUST  GAS  SPRAY  COOLER 


IMPLICIT  REAL*8(A-HtO-Z> 

COMHO»rRAtVLBtAA,lWC - - - 

COMMON  R0L,VW,GW,GC,RVtCJ,NNCA0,B(4)  ,ISTA 

COMMON  IS6T,K<K,KOUNT,KN.KKl  " 

_ COMMON  PL.BL.CL.TSA _ 

COMMON  FLN,VLA,DA,FLA 

COMMQft  GE I tGE?j^GE3jGEA »GE5»G£6iGE7j-G£8iGE9iGE10>GE II .0612,06 1  3>G61 
lA.GE15,G616kG6l7,GE18 
COMMOftJ^JLLiMA.SCFiREiAfllijmL  . 

COMMON  ROP 

_ COMMON/ LOW/.  MICK  _ 

DIMENSION  Y(99),YPI99>,A5( 10),A6< 10>,A7t 10I.ARJ  10 1 , A10< 10 > ♦ XVS I  10 
-  X )  t.X HL  !  u > ,_X HFG-t  IOI  jBARHI  10)  »QM(  10 ) iAK  t J o I  « C 0 1 1 0  ) 

Z «XHV ( 10  ) 

0 IMENS  10N.-£U_LQI*YLXLQJ  *IiLLOjj  D 11Q ) .  S  TAi  1 1 1 


103  FORMAT ( 1H1 *  *8EG  IN : STATION • » I  4 , 12X . 12 HR UN  NUMBER  ,8AR> 
_ (UO _ 1_ _ 


200 

M«M+1 

XVF-O. 0 

<i«i 
aa«  rfia 

VLB* 1 .0 

J-ujt-n 

BA*.2S 

— ■ SJ-flJ-S. 

K  1  =  { 

MICH*0  —  --  ...  -  _  - . 

READ ( 5 1 51015 , EN0*222 ) ALP 

RnanAT  1  l/Utl _ _ 

RE AO ( 5  *  102  )  NST  A 

Nsis  NSTA*1  -  ...  _  _ 

REA0(5,l0tMSTAI!)*I»l.NSl) 

READ*  9*END*222  )C V* VG*  IG»P  - -  .  --  _  .. 

PcP*14<,.0 

_ ! _ mnii.ini  i  <  ci  in.  iti.mm _ 

100 

-JjQJ _ 

REA0(5»101 MVLII >. I«1 »N$TA> 

REA015«1011  1  TS  J  1  l»AcI  «_NST  AX  .  ..  _ 

READ* 5  ? 101  KOI  I >,I*i,NSTA? 

R£A0(5*10il(B(I),l*i»*l  -  _ 

FORMAT(4El0.2> 

FORMAT  (  flPlfl.m 

102 

FORMAT*  12) 

CPL=1 .0 

GEl=STA(2)-.2500 

GE2-GE 1 ♦ . 100 

GE  3*ST A  *  3  >-. 2500 

GE4=GE3*.100 

GE5*STA(4>-.2500 

GE6aGE5*«100 _ 

GE7=STA(5>-.2*00 

GEB  =  GE  7* ♦ 100 

GE9=STA(6>-.2500 

GE 10=G£9+ .100 

GE11=STA* 7)-. 2500 

GE 1 2*GE 1 1 ♦. 1 00 

GE13*STA<8)-.25D0 

GE 1 4  =  GE l 3*. 1  DO 
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GE15*STAI9  )-.i'500 
GE16*GE1S*,100 
GE17-STAI 10>-.2*Q0 

:eib-gei7*.idq _ 

FLN-0.0 

VLA«0.0 . . . 

04*0.0 

FLA-0.0  _ 

Cl*0 .0 

_ QL*0,0 _ 

fSA-0.0 

BL*0.0  _  _ _  _ _ 

C  J  =  778 .0 

GC«  32,2  _  __ _ 

I  SET  ■  9 

_  GW=29,0 _ _ 

VW=  1 8  ,0 

RV=1  986*778.0  _  _ 

RVAP=RV/VW 

RNC-RV/GH  _  _ 

ROL-62.4 

_ TREE  *  540,0 _ 

ISTA=1 
MST  A=  t  ST  A 
N60  =  8 

C  THE  ARRAYS  USED  IN  0IFF6  ARE  NOW  SET  GP 

YIJ)=CV 
Y(2)=VG 

Y|5)=TG 

Y(4)=P 

Y(5>*FU  ID 
Y(6>=VL ( 1  ) 

Y(7)=TSI 1 > 

Y!8>=  0(1) 

- wncao  *  v r? )  *ii.o-vnYrmrwiYrmgv*m-v'nn/GW»mi7vur  r 

X  =  0.0 

WRITE(6,103)ISTA,  ALP 
3  CONTINUE 

OX®. 000 1 
KOUNT  =0 

- rrmTA  ,w.  i)j»s - 

00  1  11=1 ,9999999 

I F ( I  ST  A  .EO.  1 >KK I =8 

1  F  (  I  S T  A  .60.  2  )KK  1  =  16  __ 

IFIISTA  .60.  3  >KK I =24 

IFUSTA  .60.  4  )  KK  I  =  32  __ 

I F ( I STA  .60.  5 >KKI«40 

IFIISTA  .60,  6)KK1»48 _ 

IFIISTA  .60,  7  >  K <  I  =56 
IFIISTA  .60.  8 )KK I =64 
IFIISTA  .60.  9 )K< I =72 

I F ( X  . G6 .  S  T  A ( 2  )  ,ANO.  X  ,LE.  G63)KK|=8 
I F I  X  .G6.  G63  .ANO.  X  .LE.  STA I  3  )  KK  I -16 
I F  I  X  ,GF.  ST  A  I  3  )  .ANO.  X  .LE.  GE5>KKI=16 
I E ( x  .GF.  t65  .AND.  r  .LE.  ST  A  I  4  )  >  KK  I  =24 
I  F I  X  .06.  S  T  A ( *  )  .ANO.  X  .LE.  GE7)KKI=24 
:F(X.GE.STA(5).ANO.X.L6.G69)K((I  =  32 
I F I  X  .06.  S  T  A  I  6 )  . ANO .  X  ,LE.  GEll)KK[*40 

I F  ( x  .G6 .  STAI7)  .AND.  X  ,'.6,  GE13>KKI*48 

!  F  |  X  ,GE.  S 1  A  I  8  )  .AND.  X  ,L6.  GE15)KK|=56 _ 

I F  (  X  .06.  STAI9)  .AND.  X  ,'.6,  G617)KK!»64 
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I  F  I  I  ST  A  .60.  2  .AMO. 
1P(X  .06.  STAI2)  .ANO 
I  P (1  ST  A  .60.  3  .ANO. 
IFIX  .06.  STA ( 3 )  .ANO 
I P (  I  ST  A  .60.  A  .AND. 


X  .66 *  STA(2> (MSTA-3 
.  X  •  LE.  GE3  >  MST  A»  3 
X  .G6 .  STAI3) >MSTA«3 
.  X  .1.6.  GE5>MSTA»3 
X  .06.  ST  A  I A  > )HST  A«3 


1 F ( I  ST  A  .60. 

5 

.AMO.  X 

.GE. 

STAI5) >MSTA«3 

I F ( I  ST  A  .60. 

6 

.AMO.  X 

.GE. 

STA(6) >M$TA*3 

IF  1 1  ST  A  .60. 

7 

.AMO.  X 

.GE. 

STAC7 » )MSTA*3 

IFUSTA  .60. 

8 

.ANO.  X 

.GE. 

STA 1 8  > )MSTA»3 

IF ( 1 STA  .60. 

9 

.AMO.  X 

.GE. 

STA 19  > ) MS  T A«3 

LMX  .GE.-SFJ 

i  » 

ANO.  X  . 

Li. 

GE4)MSTA-1 

IFIX  .GE.  065  .AND.  X 
!F ( X  , GE.  G67  .AND.  X 
IP(X  .GE.  G69  .ANO.  X 
I F  (  X  . G6»  0611  .AMO. 

I F ( X  .GE.  GE 1 3  .AMO. 

I F ( X  .GE.  GE15  .ANO. 


•  LE.  G66 )NSTA» 1 
«L£.  G68)«STA=l 
.LE.  GE 10  >MST A»  1 
X  .LE.  GE12 )MSTA*1 
X  .LE.  GE14>MSTA=1 
X  .LE.  GE 1 6 ) MS  T  A«  1 


I F  ( X 

.GE. 

GE  1 7  .ANO. 

X  .LE 

.  GE  18 )MST  A* 1 

I F  ( X 

.GE. 

GE2  .ANO. 

X  .LE. 

STA(?> )MSTA=2 

I  F  (  X 

.GE. 

G64  .ANO. 

X  .LE . 

STA ( 3> )MSTA=2 

I  F  ( X 

.GE. 

G66  .ANO. 

X  .LE. 

STAI4) )MSTA*2 

I F  I  X 

.GE. 

GE8  .AMO. 

X  .LE. 

ST  A ( 5 ) )MSTAe2 

I  E  ( X 

•  GE. 

GE  10  .AND. 

X  .  LE 

.  ST  A ( 6 ) ) MS  T  A* 

I F  (  X 

.GE. 

GE  12 

.ANO.  X 

.LE. 

ST  A ( 7 ) >MSTA  =  2 

I F  ( X 

.GE. 

GE  14 

.AMO.  X 

.LE. 

ST  A ( 8  > ) MS  T  A=2 

I  F  ( X 

.GE. 

GE16 

.AND.  X 

•  LE  • 

STA(9) >MSTA=2 

IFIX 

.GE. 

GE18)MSTA=2 

CALL  OIFF6<X.Y,Yf>,OX,tKK,NEO,Kl  ,MSTA> 


KI*1 

I F  ( X 

.GE. 

GE1 

.AND. 

X 

•  LE.  GE l*OX  >G0 

TO 

909 

I E  (  X 

•  GE. 

GE3 

.ANO. 

X 

.LE.  G63+OX )G0 

TO 

666 

I  F  ( X 

.GE. 

GE5 

.ANO. 

X 

•LE.  GE5*0x )G0 

TO 

666 

IFIX 

.GE. 

GE7 

.ANO. 

X 

.LE.  G67+OX )G0 

TO 

666 

IFIX 

.GE. 

GE9 

.ANO. 

X 

.LE.  G69+0X )G0 

TO 

666 

IFtX  ,G6 .  G611.AN0.  X  .LE.  GEll+DXIGO  TO  666 


I F ( X  . GE.  GE 1 3  .ANO.  X  ,L6.  GE13*DX)G0  TO  666 _ 

IF(X  .GE .  GE IS  .AMO.  X  .LE.  GE15*0X)G0  TO  666 

1F<X  .GE.  GE17  .AMD.  X  .LE.  GE17*0X)G0  TO  666 

GO  TO  669 

666  FLN«mLA)fY(LA-<i).Y(LA-8>  > 

Vl.A=(Y(LA*l  >*Y  (LA  )*Y(  LA-3  )*Y(LA-P)*Y(  LA-7)  *Y(  LA-8)  >/FLN 
0AslY(LA+3)«YILA>*Y(LA-l)«Y(LA-4 ) ♦ Y ( L A-S ) » Y ( L A-B ) )/FLM 
TSA=( YILA+2 )*Y IL A  >*Y I L  A-2 ) * Y ( L A— 4 ) ♦ Y I L A-6 ) • Y I L  A— 8  ) l/FLN 
GO  TO  669 

909  FLN  =  FLN«-Y<  J-3) 

VLA  =  VLA+( Y( J-3  >*Y( J-2  )/FLN> 

DA=OA*( Y( J-3)»Y( J l/FLM) 

TSA«TSA*Y( J-l ) / 1 S  T  A _ 

669  CONTINUE 


IFIX 

.GE.  GE 1 

.ANO. 

X 

.LE. 

GEl+OX )G0 

TO 

8000 

IFIX 

.GE.  GE 3 

.AMO. 

X 

•  LE. 

GE3*DX >G0 

TO 

8000 

IFIX 

.GE.  GE5 

.ANO. 

X 

.LE. 

GE  S+  OX  >G0 

TO 

8000 

IFIX 

.GE.  GE7 

.ANO. 

X 

.LE. 

GE7+0X (GO 

TO 

8000 

IFIX 

.GE.  GE9 

.ANO. 

X 

.LE. 

G69+OX (GO 

TO 

8000 

I F  ( X  .GE.  GEll.AND.  X  .LE.  GE11*0X(G0  TO  8000 


I F  <  X  .GE.  GE 1 3  .ANO.  X  .LE.  GE13*0X)G0  TO  8000 

I F ( X  .GE .  GE 1 5  .AMO.  X  .LE.  GE15-*DX)G0  TO  8000 

1 F ( X  .GF.  GE 1 7  .AMO.  X  .LE.  GE 1 7+DX  >G0  TO  8000 

I P ( X  .GE.  GE2  .AMD.  X  .LE.  GE2*DX)G0  TO  8001 

m X  .GE.  GE4  .AMO.  X  .LE.  GE4*0X)G0  TO  8001 

IF  ix  .GE.  GE6  .AMD.  x  .LE.  r,C6  +  0X>G0  TO  8001 
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I  F  I  X  .GE.  G68  .ANO.  X  ,L6  •_  G68  +  OX  >00  TO  8001 
I  (=  ( X  .GE.  G610  .AND.  X  ,u6.  GeiO*DX>GO  TO  3001 

1FIX  .GF.  GF12  • ANO. X  .LE.  GE12+DX1G0  TO  8001 

I F ( X  .GF.  G614  .AND.  X  .L6.  GE14*0X)G0  TO  8001 

_ I F ( X  .GF.  GF 1 6  .ANO.  X  .LE .  GE16*DX )G0  TO  8001 

l  F  I  X  .GE.  GE 13  .ANO.  X  .LE.  G618*DX)GO  TO  8001 

99  IFIX.GE.STAI ISTA+l > >G0  T0_2 _ _  _ 

GO  TO  1 

8000  CONTINUE  _ _ 

I F ( l ST  A  .60.  1 )LA  =  5 

_ 1 F ( I S T A  «E0.  2 )L A=I 3 _ 

IFIISTA  .EO.  3>LA=2l 

IFUSTA  .EO.  4)LA=29  _ _ 

IFIISTA  .EO.  5ILA=37 

IFIISTA  .EO.  6  )L A=45  _ __ _ 

IFIISTA  .EO.  7 )LA=53 
IFIISTA  .60.  8>LA=61 
IF!  I  St  A  .60.  9)LA=69 
Y I L  A+ 1 )=VL A 
FLA=AA*FLN 
Y I C  A  )  =  FL A 

RG  =  Y ( 1 ) *RVAP+ I l.O-Y I  1 ) )*RNC 
R0G=Y(4)/lkG*Y<3»  > 

- rrT-  7Mzrm - 

YILA  +  2  )  =  TSA 
V(LA+3)=0A 
Kl«0 
GO  TO  1 

8001  CONTINUE 

fTl t  STA  .EO.  I  >  JJ  =  8 
IFIISTA  .60.  2  )  J  J  =  16 
I F ( (5Ta  .60.  3  I J J  =  24 
IFIISTA  .EO.  4  1 J  J  =  32 
IFUSTA  .60.  5  >  J  J  =  40 
IFIISTA  .60.  6  >  J  J  =  48 
IF!  I STA  .60.  71 J  1  =  56 
IFIISTA  .60.  8  1  J J  =  64 
IFIISTA  .60.  9»JJ=f2 
P I  6  =3.1416 
00  88  1=1,2 
J  =  JJ  +  4*< 1 -1  ) 

GO  TO  I  20 . 2 1 ) ,1 

20  CONTINUE 
FL A= AA*FLN 

F L 8=  I  1  .O-AA  )  »FLN 

RG  =  Y I  1 )*RVAP*( 1.0-YI 1 > >*RNC 

R0G  =  Y 1 4 ) / I RG*Y I  3  )  ) 

SIG=. 004790 

_  0B= 1 4 1 8. 6»  S IG)/IROG«(Y(2)-VL8>«»2) _ 

GO  TO  22 

21  CONTINUE 

Y ( J-2  )  =  VLB 
Y I J-3  >  =  FLB 
Y( J-l  )  =  T S A 

_ Y  I  J )  =  OB _ 

I KK  =  I KK  *8 
NEO  =  NEO  +  ‘‘ 

K  [=0 

22  CONTINUE 
88  CONTINUE 
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1  CONTINUE _ 

JJJ  =  JJJ*l 
2_  _  CONTINUE 

WR1T6I6, 12139) 

12139  FORMAT  ItHl  ) 

C  we  WIU.  NOW  PLOT  THE  STATION  JUST  FINISHED 

_ I_S_T  A«  I  ST  A*1 _ 

IF  (  I  ST  A  .60.  1)U'5 
I F ( I  ST  A  .EO.  2  )l As  1  3 
IF ( I  ST  A  .60.  3>LA=21 
IFIISTA  .60.  4>LA=29 
IFIISTA  .60.  5>LA=37 

_  IFIISTA  .60.  6)LA=4S _ 

IFIISTA  .60.  7  )  L A=  5  3 
IFIISTA  .60.  8>LA=61 
I F  I  ISTA  .60.  9)LA=69 
IFIISTA  .GT.  NSTA)GO  TO  5 
I  SET  s  9 

_ WRIT6(6,103)1STA,ALP _ 

C  SET  UP  ARRAYS  FOR  DI6FE 

YINE0»1  )  =  FLI  ISTA) 

VINEO+2 >=VL I ISTA) 

YIN60*3)=TSI ISTA) 

Y(NE0+4)=0I ISTA) 

_ N60  =NEQ*4 _ _ _ 

KflUNT  =  0 

GO  TO  3 _ 

5  GO  TO  200 
222  STOP 
ENO 


SUBROUTINE  OIFF£|X,Y,YP,DX, IKK.NEO.KI ,MSTA  ) 

IMPLICIT  REAL«8I A-H.O-Z I 

COMMON  BA.VL8,  AA.RNC _ ; _ 

COMMON  RGL  .VW.GW.GC.RV .C  J. WNC AO. 8(4)  .ISTA 

COMMON  IS6T.KKK  .KOUNT .KN.KK I 
COMMON  OL.BL.CU.TSA 
common  FLN.VLA.OA.FLA 

COMMON  GE1 .G62,G63.G64,G65,GE6,GE7,G68,G69,GE10.GEH.GE12,G613.G61 

14.GE15,GF16,Gf  17.GEI8 _ _ 

COMMON  PVTS.OAB.SCF.RE.ABN.FNU 
COMMON  ROF 

DIMENSION  YI99)  ,  YPI99) ,Z  <99  > ,ZP199>,  ZNI 99) 

DY=0.0 

CALL  YFUNCIX.Y.YP.KI  ,  l , MS T A , I KK , OX , OY ) 

120  CONTINUE _ _ 

00  1  1  =  1.  N60 

Zf I  )  =  Y(  I  ) +0  X  *  YP  I  I  ) 

1  CONTINUE 
X=X+DX 

0x2= .50*0*OX 

00  5  J  =  2.999 _ _ _ 

OY  =  0X 

666  CALL  YFUN'C  I  X,  Z  ,  ZP.K  I  ,  J.MSTA  ,  IKK  .OX.OY  ) 

K  =0 

00  it 0  1  =  1  ,NEO 

ZNII  )  =  Y ( I  )  *0X2* I YP  I  I  )  »Z  P  I  I  >) 

IFIOARS  I  ZNI  I  )-Z  I  I  )  )-l  .0-05*0ABS(ZNI  I  )  )  )4,4,  3 _ _ 

3  iT=  1 


59 


AEOC-Tfl.72-89 


4  KK  =  J 

l C I )«2N( I > 

40  CONTINUE 
IE(K>5,6,5 

_  5  CONTINUE _ 

WRITE  (  6.99  K2N<  I),I=1,N601 
WK I TE < 6.90  )KK 
90  FORMAT  I  15) 

99  FORMAT (4620.10) 

STOP 

6  00  7  1=1. NEO _ _ 

TTTT )  =  Z  (  I  ) 

1210  CONTINUE 

I F ( J  ,GE.  3  „ANO.  J  .LE.  5)G0  TO  1212 
IF ( J  .Li .  3)00  TO  2020 
0X=.5*DX 
GO  TO  1212 

2770'TSTiT70^TTx - 

I F ( OX  .GT.  .01 ) DX= . 01 

1212  CONTINUE 
RETURN 
6N0 


SUBROUTINE  YFUNCIX.Y.YP.KI .K.MSTA, IKK.DX.DY) 

- ,HPL  rcTTBFTrr*!  i  A-H-;n-z  > - 

COMMON  BA. VLB, AA.RNC 

Common  R0L.vw.6w,ftC.RV,Cj.wN£A6,8<4l , , iSta 

COMMON  ISET.XKK.KOUNT.KN.KKI 
COMMON  ol.bl.cl.tsa 
COMMON  FLN.VLA.OA.FLA 

COMMON  f.El , GE2,G63,GE4,GES.GE6,GE7.C>6fl,G69,GE10»GEH»GE12.G613,GEl 
14.GE15.GF16.GE17.GE18 

Common  pvts ,oab, SCf.re , a'bn.Enu  ! 

COMMON  ROF 

DIMENSION  Y(99 ) ,YP(99) , A5( 10) . A6( 10) , A7I 10) . A8( 10) , A 10 1  10) , XVSI 10 
1  )  .XML  (  10)  ,XHFG<  10)  ,BARH<  10)  ,0M(  10).XK|,10)  ,C0(  10) 

2,x  hV (10) 

RVAP=RV/VW 

WJ^'TMTTOTE - 

I F (  ISTA  .EQ.  I  )  J  =  8 
I F (  ISTA  .EQ.  2  >  J-  16 
I F ( ISTA  .EQ.  3  >  J - 24 
I F ( ISTA  .EQ.  4 ) J-32 
IF (  I STA  .EQ.  5  1  J  =  40 
IF (  1 STA  .50.  6  >  J  =  48 
I F I ISTA  .EQ.  7) J=56 
I F (  | STA  .EQ.  8  )  J  =  64 
IF! ISTA  .EO.  9)J=7 2 
SIG=. 00479 

_ PIE  -3.1416 _ _ 

SUM  1=0. 

SUM2=0.  i 

SUM  3  =0  . 

S  UM4  =  0  # 

SUMS  =0 . 


I  F  {  X 

•  GE. 

197cm 

■fum 

JL 

•  LE. 

9939191 

171 

II'MI 

IF  (  X 

•  GE  . 

GE2-0X 

•  AND. 

X 

.  L  6 . 

GE2  )G0 

TO 

8001 

I  F  (  X 

•  GE. 

GE3-0X 

.ANO. 

X 

•  LE. 

GE  3 )G0 

TO 

8001 

I  F  (  X 

.GE. 

GE4-DX 

•  M  *»  # 

K 

.LE. 

GE4  ) GO 

TO 

8001 
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IF<x  .GE.  GE5-DX  . AND*  X  .LE.  GE5IG0  TO  BOOl 
IFIX  .GE.  G66-DX  .ANO.  X  it E.  GE6IG0  TO  8001 


IF(x  .GE. 

.LE. 

irTlJlfTiWTr 

8001 

1 P I  X  .GE. 

G68-0X  . 

AND.  X 

.16. 

GE8IG0  TO 

8001 

I F  I  X  .  Gc  • 

GE9-0X  . 

ANO.  X 

.LE. 

GE91G0  TO 

8001 

I F  I  X 

G610-0X. 

ANO.X.LE.GElOlGO  TO  8001 

I F ( X  .  E. 

GFll-D* 

.ANO.  X 

.LE. 

GE111G0 

TO  8001 

I F ( X  „£. 

GE 1 2-OX 

.ANO.  X 

•  LE. 

GE  1 2  >  GO 

TO  8001 

|F(x  .GE. 

GE13-0X 

.ANO.  X 

.LE. 

GE13IG0 

TO  BOOl 

I F ( X  .GE. 

GE14-0X 

.AND.  X 

.LE. 

GE 14 ) GO 

TO  8001 

IF ix  .GE. 

GE 1  5-OX 

.ANO.  X 

.LE. 

GE I 5  IGO 

to  eooi 

I F ( X  .GE. 

GE 16-OX 

.ANO.  X 

.LE. 

GE 16 ) GO 

TO  BOOl 

I F ( x  .GE. 

GF 1  7-OX 

.ANO.  X 

.LE. 

GE17 IGO 

TO  8001 

I F l X  .GE. 

GE 1 8-OX 

.ANO.  X 

.LE. 

GE 18 IGO 

TO  8001 

9998  CONTINUE 


8001  CONTINUE 

RG*Y< 1 l*RVAP* ( 1 .0-Y< l ) l*RNC 

8002  IFIKI  .60.  0 )G0  TO  8003 

CALL  ROGEE (ROG.Yl l I ,YI21 ,X I  , 

_ GO  TO  6009 _ 

8003  R0G=Y<41/CRG*Y(3> > 

8009  CONTINUE _ _ 

00  l  1*1 tHSTA 
J«KKI>4*( l-l » 

AS  (  1  >  *  P  1 6*R0L* Y ( J  1**2  /2.0 

CALL  SUB1  CYll>,YI3l,YI4l,Y<J-l),XVSU  I.XV.lfKin.YI  JI.DMI  I  I.CDI  I  >.B 
IARHI  I  ),Y|J-2I,Y<21> 

_ |F(XVf  .GT.  1 .0  I  RETURN _ 

VP(J>=  mm*  PIE  *Y(J)**2  »I  XV-XVSI  I  >  )*VW/(  A5I  I  >*Y(  J-2  )*(  1-XVSC  I  ) 
111 

A6 (  I  ) =  VP ( J  > 

YP< j-2  )  =  < <P [E/8.0 l*ROG*< Y<  Jl**2  l*CD( I  I *DABS IY( 2 )- Y < J-2 I )»( Y( 21- 
1 Y  (  J-2  I  )-(  Y(2>- Y( j-21 l*Y< J-2)*AS<  I  >*A6I  m/IOMt  I  )*Y<  J-2II 

A7( I )=YP( J-2  ) _ 

YPU-31*  A5(I l*A6( I >  *Y I J-3  » /DM  I  I > 

A8( I )*YP( J-3) 

SUM1  =  SUM1*A5( I )*A6( I l*Y( J-3I/DMI  I  ) 

CALL  HL! ( Y( J-l ) ,XML ( I >  »  XHFG ( I  > , XHV I  I  ) > 

YPIJ-1)  *  B  ARH (  l >  *  PIE  *  Y(J>*»2  *< Y< 31-YI J-l ) I /  I  DM!  I  ) * Y < J- 2  )> - 

CXHL (  I  ) «AS I  I l*A6(  I  ) /QM( I  > ♦ A 5 ( I  I  *  A  6  (  I  I  ♦  X  H  V I  I  > / QM 111 _ 

A I  0 (  I  ) =  YPI J-l ) 

SUM2  =  SUM2  +  YU-’  •  i  i  »♦  Y  (  J-3  l*A7  (  I  ) 

SIJM3  =  SUN  3  *A8<  i  >*<XHLU  l+YI  J-2)**2  / ( 2 . 0*GC*C J I  I 
SUM4  =  SIIM4  * Y (  J-3  i *  Y(  J-2  l«47l  I  l/(GC*CJ) 

SUMS  =  SUMS  ♦  Y I J-3 ) *A1 0( I > 

1  CONTINUE _ _ 

YPIll  =  -SUMl*( 1  — Y I  1  )  1**2 
A9=»YP ( l 1 

CALL  AAA3I X , Y( 1 ) , Y < 2  I » A3, AAO.DAAO.ROG) 

CALL  OERT (  RUG , Y ( 1  I , Y ( 2 ) , Y ( 3 ) ,  Y ( 4 > , A9. SUM2 . SUM3, $UM4, SUMS . AO , A  3 . A  1 
C l, A l 2, A l 3. A  1 4, R, CP V, CP A, A 2, A4, A AO. A) 

_ YP<3>  *IA11«A14-A12»A13)/(AH»(CPV»YI1)/(1-Y(1)  I  *C P A > - A  1 2* AO»RQG» R 

Cl 

A1S=  YPI 3 ) 

YP (2)  =  I A13-AO«ROG«K*A15 I /All 

A l 6  =  YP ( 2  I 

CALL  06 RP ( RCG.YI 11,YI2I,Y(3I,R,A2,A3,A4,A9,A15,A16,0P,A) 

_ YP  I  4  )  =  OP _  _ _ 

IMK.GT.  1)  60  TO  3 
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IFIX  .60.  0.0)G0  TO  9999 
1 F ( X  .IT.  64  >G0  TO  3 
BA-8A+.25 
9999  WRITE<6,4)X 

_ WRUEI6,fi8888)R0F 

"SW«rll  )/(  l.O-V  (  1  )  » 

4  FORMAT I lHO.  • X« • ,616.10) 

WRIT6I6.S*  4A0, WH 

"5 - FORMAT!  1h6‘,  'i/iO*'  ,Urr.“675)f P ffV/wNC"-  ‘  ,614.61 - 

WRtTE(6,6>Y(l>»Y<2),vm»Y<4> 

I F ( 1  $ T  A  .FQ.  I >  J  J  =  8  • 

IFUSTA  .00.  2  > J J*  1 6 
IFUSTA  .60.  3  >  J J  =  24 

_ IFUSTA  .60.  4)JJ»32 _ 

lFdstA  .to.  5)jj=4"o 

IFUSTA  .00.  6  )  J  J  =  48 
|F(  I  St A  .60.  7)  JJ  =  56 
IFIISTA  .00.  8  )  J  J  =  64 
IF! ISTA  .60.  9)  JJ  =  72 

6  FORMAT (lHO. *CV=* ,G20.9,SX,*VGs*,G20.9,5X, • TG» • , G20, 9, 5X , 'P«< .G20.9 

l  > 

_ IFCX  .06.  G62  .AND.  X  ,L6 ♦  G63-OX)GO  TO  88 _ 

1 F ( x  .06.  GE4  .AND.  X  .16.  G65-0X)66"T0  88 

I F ( X  . G6 .  GF6  .AND.  X  .16.  G67-OX>GO  TO  88 

I F ( X  .06.  G68  .AND.  X  .L6 .  G69-0X  >G0  TO  88 

IF(x  .06.  G610  .ANO.  X  ,L6.  G6U-0XJG0  TO  88 
I F ( X  . G6 .  G612  .ANO.  X  ,L6.  G613-OXIGO  TO  88 

_ I F ( x  ,oe.  6614  .ANO.  X  . L  6 .  G6I5-0X)G0  TO  88 _ 

1 6 ( X  . G6 .  0616  .ANO.  X  ,L6.  G6 17-OX  >G0  TO  88 
IFIX  .06.  G6 18-OX >G0  TO  88 

on  io  1*1,1 

J=JJ*4* ( |-11 

WRITe<6,fl>  y( j-3 > ,V( J-2 ) ,Y( J-l  ),Y f J> 

_ WRITE <6,91  YP(  J-31.YPIJ-2)  ,  VP  I  J-l  >  ,  VPI  J  ) _ • 

10  CONTINUE 

GO  TO  1001 
88  00  988  1  =  1  ,2 

J= I KK  +  4* (  1-1  1 

WR I T6 Ife.R )Y( J-3) ,Y< j-2) ,YI J-l >,V<  J) 

WR|T6(fe.9)YP(J-3>.YP(J-2).YP(J-l).YPCJ) _ 


988  CONTINUE 

1  6  (  X 

.06.  OF  1 ♦ . 2600-Ox 

.AND. 

X 

.16.  G63-OX 100 

TO 

1000 

IF(X 

.06.  06  3*- .  2  SOQ-DX 

•  ANO. 

X 

,L6.  G65-DX )G0 

TO 

1000 

IFIX 

.06.  OF  5* . 2  500-DX 

.ANO. 

X 

.16 .  GE7-0X1G0 

TO 

1000 

IF1X.G6.  GF7+.2500-0X.AN0.  X  .UE.  G69-0X1G0  TO  1000 


IFix.06.  GF9* . 2  5000-DX  .ANO.  X  . 1 6 .GE 1 1 -OX >00  TO  1000 


IFIX 

.06.  G6 lt+. 25O0-OX 

.AND. 

X 

.IE. 

GE 1  3-OX  )G0 

TO 

1000 

IFIX 

.06.  G61 3+.2500-0X 

.AND. 

X 

.16. 

G615-DXIGO 

TO 

1000 

IFIX 

•OE.  06 15* .2SD0-DX 

.AND. 

X 

.IE. 

GE17-DX  >G0 

TO 

1000 

GO  TO  1002 

iooo  Continue 

_ I  F ( X  .06.  061 -OX * .2800  >  J J  =  8 

I  F ( x  .06.  G63-OX+.25DO) JJ=16 
I  F  (  X  .GE.  GF8-0X  +  .2500) JJ  =  24 
I  F ( X  .06.  OF7  +  .2SOO-OX ) J J=32 
'FIX  .06.  GF9+.25no-f)X  >  J  J  =  40 
I F I  X  .06.  OF  1 l  +  . 7800-OX  1  JJ  =  4R 
_ I  FIX  .06.  GF13*.250Q-0X)JJ=S6 

7  =  i x  .oe.  06  is < . 2'  'O-ox ; j j=64 

IFIX  .OF.  OF  1 8— OX  >00  TO  100? 
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1010  J-JJ+8 

WRITE (6, 8)  Y!J-3),Y!J-2>,Y< J-l>,Vf J> 

WRITE  16,9)  YPIJ-3) ,YP( J-2 > , YP ! J-l ) , YP t J) 

199  CONTINUE _ 

1001  CONTINUE 
1002  CONTINUE 

WR1TE<6,7)YP( 1 > ,YP(2) ,YP<  3) ,YP(4) 

8  FORMAT! IHO, 'FL»* ,G20.9,5X, 'VL** ,G20.9, 5X, 'T$«' ,G20.9,5X, 'O*' ,G20.9 

1  > 

9  FORMAT! 1H0, ♦FLP«« .G16.6.5X, 'VLP* 1 , G16. 6, 5X , » TSP- » , G 1 6. 6 , SX, » OP* « , G 
116.6) 

7  F0RMAT11H0,,CVP3,»G16.6,5X,»VGP**,G16.6,5X»>TGP»>,G16.6.SX, »PP*»,G 

116.6) 

3  CONTINUE 

2  FORMAT  1 1H  , 7E 16 .6  > 

RETURN 

ENO 


_ SUBROUTINE  SU81!  C Y , TG , P , T S . X V S . X V , XK , 0. DM, CO. BN, VL . VG  ) _ 

IMPLICIT  REAL‘8! A-H.0-2 > 

COMMON  BA ,  VLB . A A , RNC 

COMMON  R0L,VW,GW,GC,RV.CJ.WNCA0.B<4>  ,ISTA 

COMMON  I$FT,KKK,Kf)UNT,KN.KK| 

COMMON  OL.B'-.CL.TSA 

_ COMMON  eln,  VLA.OA.FLA _ _ 

COMMON  GE 1, GE 2, GE 3, GE A. GE5. GE6, GE 7 , GE8, GE9, GE 10, GE1 l, GE1 2. G61 3, GE1 
14,GE15,GE16,GE17.GE18 
COMMON  PVTS,DA8,SCF,RE.ABN,FNU 
COMMON  ROF 
9899  CONTINUE 

_ XV=IC//VW)/(CV/VW  ♦ I  1 -C V ) /GW  ) _ 

X2=-9.06*( (-.5696*T5+  .00396-4  *TS**2  +.0927E-7  *  T  S**  3*  1  352 . 3 )  /  T  S- 
C1 .9425 ) 

X 1=  672.0/TS 

PVTS  =  2117.0  *Xl**5.19*06XP!X2> 

XVS  =  PVTS/P 

XVF  =<XVS*XV >/2.0  _ 

X Ms  XVF*  VW  -Ml-XVF)  *r,W 
TFI  =!TS*TG)/2.0 

rMV  = ( 10.6E-7  *OSORT(TFI  ))/!  1  +  1538. O/TFI I 
FMNC=  7.5E-7*OSO«T(TFI )/( 1+  216. O/TFI) 

FM  =  XVF  *  F  MV  ♦ (  1 -XVF ) *FMNC 

OAfi  =  ( ,5375/P >* ( TF I /491 .0 ) **2. 334 

ROF  =  P»  XM/I RV*  T  F  I  ) 

SC  F  s  FM/ (RDF  *0 A  8 ) 

R6  =DABS(  VG-VL  >*Rf7F*D/FM 
IFIRE  .LT.  0.0 ) GO  TO  9999 
4BN=  2.0+0.6*DS0RT (RE )  *  SCF**  0.3333 
XK  =  ARN*FM/ ( 0*SCF *XM ) 

OM  =  0**3  *3.1416  *ROt»  /o .  0 

CO  =  24.0/RF*!  1 .0*. I 5*RE**0.687  ) 

I F (  TF]  ,GT. 1700.0  .AND.  TFI  .LT. 4500.0)  GO  TO  221 
IFITFI  .LT.  400.0  .OR.  TFI  .r,T.4500>  GO  TO  222 
CPVF  =.4304  ♦•  .  167BE--.*TF  I  +  0 . 2  78  i  E- 7*  T  F  I  *  *  2 . 0 
C PNC F  =  0.2318  ♦  0.1040E-4  ♦TFI  ♦  0.7166E-8  *TFI**2 
GO  TO  102 

221  CPVF  =  .3319  ♦0.143HE-3  *TFI  -0.1312E-7  *TFI**2.0 

- CP7JC?="  07?FU  -T  O'.  J5?1E-4*TF'I-  0.  3776T-S"  *771**5 - 
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GO  TO  102 

222  WRITE(6,101> 

101  FORMAHlH  ,  •  TEMP-F  I  •  I  S  OUT  OF  RANGE*) 
WRIT6<6,122>  TG.TS 

122  FORMAT(lH0**TG“,tE20.6«  5X«  *TS«*  ,620.6) 

vq^VCgNTTNUE - 

_  WRITE(6.9990)RE.VG.VL.R0F,D.FM _ _ 

"  999?TE0RMATl6EiB.10> 

STOP 

102  CONTINUE 

FKV  a  0.432*FMV 

FKNC  «  0.257  *(  0.  I 15*  5.17*CPNCF>*  FMNC 

F  K  ■  XVF*FKV  ♦(1-XVFI*  FKNC  _ _ _ 

tW  »  xvf*(vw/xM  >»£Avf  ♦( 1-xvF )*tP^CF*Gw/xM 

PVF  a  CPF  *FM/FK 

FNU  =  2*. 60  *DSORT ( BE )  *  PVF*a  0.3333 
BH  =FNU  *FM  *CPF/(0*PVF) 

RETURN 

End 


subroutine  RoGEEtrtOG.CV.VS.STi 

IMPLICIT  RFAL»8( A-H.O-Z  ) 

COMMON  BA,VLB,AA,RNC 

common  ROL.VW.GW.GC.BV.CJ.WNCAO.BIA)  ,ISTA 

& AO  =BI1)  ♦  8  (  2  )  *  X  +B(3)*X**2.0*  B(4)*X**3.0 
ROG  =WNCAO/( < l-CVl*VG*AAO> 

RETURN 

ENO 


SUBROUTINE  HLKTS.XHL.XHFG.XHV) 

IMPLICIT  RE AL  *8  I A-H ,0-2  ) _ 

XHL  =  TS  -540.0 

XHFG  =  -.5696*TS  +  .08396-4  *TS**2  ♦  0.0927E-7  •  TS**3. 0* 1 352. 3 

XHV=XHL+XHFG 

RETURN 

END 


SUBROllT I NF  AAA3(X,CV.VC.A3.AAO.DAAO.ROG) _ 

IMPLICIT  REAL*8( A-H.O-Z  ) 

COMMON  BA, VLB.AA.RNC  _ 

COMMON  ROL .VW.GW.GC.BV ,C  J.WNC AO, 8(4)  ilSTA 

COMMON  I  SET.KKK .KOUNT.KN.KK I 

A AO  a  B(l)  ♦  8(2  )  *X  ♦B ( 3 ) *X  *  *2 . 0  ♦  B < 4 ) *  X** 3. 0 
PA AO  =  B ( 2 )  +2.0*  8 ( 3  >  *X  *3 .0*6 ( 4  >*X**2. 0 
A3  =  VG  *<1-CV>*  OA AO/  ( VG*( 1-CV)*AA0)**2 
RETURN 
END 


SUBROUTINE  PERT (  ROG .C V . VG . TG ,  P , A9 , SUM2 , SUM3,  SUM4, SUMS, AO, A  3, A  11 
1A  12  ,A13,A14,R,CPV,CPA,A2 ,A4, AAO, A ) 

IMPLICIT  RFAL»8( A-H.O-Z ) 

COMMON  BA, VLB.AA.RNC 

COMMON  ROL , VW.GW.GC .RV.CJ.WNCAO, B(4)  ,ISTA 

COMMON  ISET  ,  KKK 

IF(VG  .LT.  O.OlGO  TO  20 _ 

R=  P/ ( ROG*  T& ) 
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AO  *  GC* ( 1 -C V ) /ROG 

AA  «  WNC AO/ (VGA*2*(1-CV)  * AAO ) 

All  «  ( 1 • 0-C V )*VG-TG*R*AO*AA 
A12  =  (1+CV/(1-CV> >*VG/(GC*CJ  > 

A 1  =  GW  »CV  »VW»(1-CV) _ 

A  =  WNCAO/(VG*(1-CV)**2*AAO) 

XV«»  GW  »C V / A 1 _ _ _ 

A2  .  «  RV*VW*GW* (GW-VW )/<Al**2*(XV*VW*( 1-XV)*GW)»*2| 

A1 3  »  A0*WNCA0*TG*R*A3-A9*< VG**2*ROG*TG*A2*AO*TG*R*A  *AO)  -VG*(1-C 
CV ) **2*SUM2 

CALL  HEEV  (TG.XHV.CPV, CPA  ) 

A 1 A  »-A9*(XHV-*-VG**2.0/(2*GC*CJ)>/<  1-C V)**  2  -SUM 3“ SUM A- SUMS 

_  REJURN _  _ _ 

20  HRI T€ (6.31  )VG 
31  FORMAT (F20. 10) 

STOP 

END 


SUBROUTINE  HEEV( TG.XHV.CPV.CPA) 

I  MPL  I CJ  T  RF  AL^B  (  A-H_»0^2J _ 

COMMON  RA.VLB.AA.RNC 

COMMON  R0L,VW,GW,GC.RV,CJ.RNCA0,B<  A>  ,ISTA 

COMMON  ISET  . KKK 

I F (  TG  .GE. 1700.0  .AND.  TG  .LE.ASOO)GO  TO  222 
I F (  TG  .LT.AOO.O  .OR.  TG  .GT.A500)  GO  TO  221 

_ XH V  a  , A30A»TG  ♦.OB39E-A»TG*«2.0«-0.0927E-7»TG»»3.0-236.3161*10A2.9 

OPV  =  . A30A  ♦  .1678E-A*T(.+  .2781E-7*TG**2.0 
CPA  =  .2318*  . 10A0E-A*TG*.7166E-8  *TG**2.0 
GO  TO  100 

222  XHV  =  .3319  *TG  ♦  . 07 1 9E -3«TG**2 . 0-.OA373E-7*TG* • 3- 1 85. 38 1 ♦ 10A2.9 
CPV  =  .3319  ♦  . I A38E -3*  t  G-  0. 1 312E-7* TG**2 .0 

_ CPA  =  «  ?? 1  A  ♦0.35?1E-A«TG  -0.3776E-8»TG*»2.0 _ 

GO  TO  100 
221  WR I T  F ( 6 • 99 ( 

99  FORMAT (  1H  ,'TEMP-G  IS  OUT  OF  RANGE') 

WRITEI6.1)  TG 

l  FORMAT ( 1H0, *TG=' ,E16. 6) 

_ STOP _ 

100  CONTINUE 
RETURN 
END 


SUBROUTINE  OERPIROG.CV.VG.TG.R, A2.A3, AA, A9, A15, A16.0P, A) 
IMPLICIT  REAL*B( A-H.O-Z  ) 

COMMON  BA, VLB. AA.RNC _ 

COMMON  ROL.VW.GW.GC.RV.CJ.WNCAO.BI A)  ,!STA 

0P=  R0G*R*A1 5  ♦ (R0G*TG*A2+TG*R*A  >*A9-TG*R*AA*A16-WNCA0»A3*TG*R 

RETURN 

END 
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APPENDIX  IV 

A  LISTING  OF  A  VARIATION  OF  THE  COMPUTER  PROGRAM 
FOR  ZERO  HARDWARE  BLOCKAGE  OF  THE  DUCT 


The  following  computer  program  was  programmed  for  the  IBM  360 
computer  and  has  been  used  to  calculate  data  for  two-phase  flow  condi¬ 
tions  with  zero  blockage  of  the  ducting  and  also  for  cases  where  a  drop- 
size  distribution  was  to  be  simulated. 
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DETERMINATION  OP  EXHAUST  GAS  COOLER 
INLET  CONDITIONS 

IMPLICIT  REAL*8(A-H,0-Z) 

REAL**  T2TW (150) 

RE AL*4  AREA,RAIR,RFUEL 
REAL*4  ARTM.ARPK.ARWl 
RE AL*4  IRUN ( 2  > , WHAT t 14 ) 

- C0HWUN~7GR~P/  BKTWT  151TI  ,'AKPR  I  15U  ) ,  AHB1 1  ISO)  .  NPU - 

COMMON  /  ENTH  /  T,  TT,  R,  C0A(71»  H(6>,  C02A(7>,  M2AI71,  XN2AI7), 

1  02AI71,  H20 A ( 7  > •  WTI61.  C0BI71,  C02ftt?>,  H2RI7),  Xti2R|7), 

2  02R ( 7  )  *  H20ft<7>.  C(6).  A(6>,  P(212> 

CALL  ERR$ET(261, 256, -1,11 

HTIl>  »  2ft. 011 

- HT(?)  -  <4Torl - 

Wt I  3)  *  2.016 

wT(4>  -  2ft  .016 
WT ( 5 )  *  32.0 
wT(6>  *  18. 016 
P ( 201 >*  11.766 

- PT7irri^-TT775T3 - 

P ( 205  )  =  12.770 
P(207>*  13.298 
P ( 209  > “  13.844 
P( 2  1 1  )“  14.408 
R«  1.9A726 

z - pres  riifft  limn — Trnn - 

c 

RE AO  (5,100)  (P(  1  ),  1*32,199  ) 

RE AO  (5,100)  (Pll). 1*200,212,21 
100  FORMAT  ( ft 0 1 0 . 0 ) 

C 

Z  rEM“PERATURrcgFTF"lCI6NTS 

C 

READ  (5,101)  ( COA (11,1=1,7),  (COB( I >,1*1,7),  ( C02AI I  ) , I  *  1 , 7  ) , 

1  ( C  02  B 1  I  >.1=1,7),  ( H2  A(I  1,1*1,71,  <H28I I  1,1*1,71,  ( XN2 A (  I), 1*1, 7), 

2  ( XN2B i I ) , I  =  1 , 7 ) ,  ( 02A ( I  1,1*1,71,  (02R( I  1,1*1,71,  |H20A(  I  ) , 1*1,7), 

3  ( H20B ( I  1, 1=1,7) 

TOT  FORMAT  (5016. 7/2016.  f) 

INPUT  CONDITIONS 
907  FORMAT  ( 20A4 ) 

MASS  FRACTIONS 

Z 

105  FORMAT  (6012.5) 

l  READ  ( 1 7  , END  =  999  )  C , P 1 , T 1 , TW 1 , VO  1 , VWl , TGuE SS , I RUN, WHAT 
READ  (17)  AREA, RAIR.RF DEL 
TGUESS=TW1*.01*(TI-TW1) 

T2MAx=  Twi-.Ol 

- TjBIHirnJTTTffT - 

IF  ( P 1 . E 0.0 .AND. Tl.EO.O. 1  GO  TD  999 

- WPTT  =75 - 

C  __ 

EM  =  0. 

00  6  I  =  1,6 

EM  =  EM  +  C (  I  )/  wT (  I  ) 

7r  crroT  tottf 

EM  =  l  . 0  /  F  M 
Cvl  =  C(6> 
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CNC1  «  1.0  -  CV1 

PVl  *<P1  *  EM  *  CVl)/  lfl. 

PNC  1  *  PI  -  PVl 

- FHNC~TT>T  g~~rr  CTTCT7  7  PTTC1 - : - 

T»Tt/  l.fl 

CALL  ENTHAL 

SUM  »  0.0 

00  7  I  *  1,  5 

SOM  a  SUM  ♦  cil)  *  H(l) 

- rcwriTiug - 

HNCl  =  SUM  /  CNC 1 

T - 1 — ITER  An;  FTTH  PTR  P  POOTC  ~  A  S'  THE  ON  CITCTN  ~0  F  "ETlTffAL  FT - 

C 

00TM«RA1R+RFUEL 

CONSTNs  I  (  1545.  *00  TM  )  /  |  28.85*144.  *P1*AREA)  )**2 
C0NSTriBCONSTN/(2.*778.*32.2) 

TGSI*T1/1 • R 

- ro  s  2^  i  t  i  -  a.~s  rr  t )  7  r h - 

T  =  TGS  1 
CALL  ENTHAL 

CALL  SOMIT  ( C , H, CNC l , ENT  1 ) 

T«T*1  .8 

FUNC 1=ENT1+C0NSTN*T*T -MNC 1 

- r,TGS;, - — 

call  ENTHAL 

call  SUMTT  ( C  »H,CNC 1 *6N 1 2  > 

T=T*l.fi 

FUNC2bENT2»CONSTN*T*T-MNC1 
AO  1  TGIF« I TG$l*TGS2>/2, 

- rR-foansrTTc.s?-TG3n/tGS2r;tr.o;io-o8)-co  to  «io - - 

T  =  TGI F 

- C  ALL  '  ENTWtP. - 

CALL  SUM l T (C*H,CNC1.ENT3) 

T=T*1.P  “ 

FUNCIsf NT3+C0NSTNBT*T-HNC1 
T S T  A  * FUNC  l  *PONC 3 
TSTP=FUNC2*FUNC3 

- IF  |  1ST  A  I  •Af7A,-A0y.~A02  — " 

402  IF  ( T $ T3  )  405.405,403 

403  WRITE  (6,411  ) 

411  FORMA T ( ' 0  NO  ROOT  FOR  INITIAL  ENTHALPH*) 

r,n  to  l 

404  FUNC2  =  FllNC3 

- IT.S2  sTGTF - 

G(1  TO  401 

405  FUNC 1 s FUNC  3 
T  G  S  1  =  T  G I  F 
GO  TO  401 

410  ENTNU=  ENT  3 

- (T=7  0T>  *TTP 

C  END  OF  ITERATION 

VGl  =  ( 00TM*1 545 ,*T 1 ) / ( P 1* ARE A* 144, ) 

VGl»VGl/2«,85 

WRITE  (6,413)  ENTNU , T1 . VGl 

WR I TF I 9)C ( 6> . VG1 ,T1 ,P1 

- 4~1T  FggMAT  ~(  Tg  ENThAl*3  V  *  ‘.Erf'.'i',*  GAS  TEMPER  ATuftE  «  '.El2.4, 

*  •  GAS  VELOCITY  =  • ,E 12,4 1 

Go  TO  1 
999  COM  INUE 
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ENOF iLf  9 

— rewruo  9 

RETURN 

€NB* 


- n^r  TWTMIL - : - 

IMPLICIT  REAL*8<A-H,0-2> 

COMMON  /  ENTh  /  T,  TT,  R,  C0M7>»  H(6),  C02AI7).  H2AI7),  XN2A(7), 

1  02AI7),  H20AI7),  WT<6>,  C0BI7),  C02B(7),  H2R(7>,  XN2B(7), 

2  028(7),  H208 ( 7  > »  C<6>,  i(6>,  PI212) 

3  TT=R*T*1.8 

— rr  n'.irvTrnoT'  on  tp  io - 

00  4  J«l,6 

-vxn  rscuxnn 

CALL  HTRT  (  r  ,4,H( 1 ) ) 

Htl  )=H(  l  )*TT/WT(  1 1  +284R8 .3  M.8/WT(1) 

00  5  J<=1»6! 

5  A  ( J  )  =  C02AI J  > 

CALL  HTRT  ( T  t  A  t  H ( 2  > > 

- W(2  )-=R7 2T*7T7VU2T“  96290 * T . R7WT !  2  » - 

00  b  J  =  1 . 6 

6  A ( J ) -h2 A ( J  ) 

CALL  HTRT  $  T  t A,h<  3  > ) 

H(3>=H(3)*TT/wT 13)  +2023.8  *1.8/wT<3> 

DO  7  J*l,6 

~7"  rr  jtxxnta  n  r - 

CALL  HTRT  ( T,A,H(4) )  1 

HI4>=H(4)+TT/wT(4)  +2072.3  *1.8/wTI4) 

00  8  J  =  1.6 

8  A ( J ) =02  A ( J  > 

CALL  HTRT  ( T,A,H(3> ) 

- RT  3  V*HI  6  )  *  TT/WT  (bT7  2TT7V.7-'*l.H/WI  IS) - 

00  9  J  =  l.6 

9  A { J ) =H20A ( J ) 

CALL  HTRT  (  T  ,  A  ,H<  6 )  ) 

Ht6)=H(6)*TT/WT{6)  +60164.7  *1.8/WT(6) 

HV2  =h ( 6  ) 

- GTT_TTrT7 - 

10  no  u  j= i « 6  j 

11  A  I J l =COB ( J  I  j 

CALL  HTRT  I T l ) ) 

HIl)*H(l)*TT/WTCl>  +28488.3  *l.B/WT(l)  3 

00  12  J=l,6  j 

TT  A”(  J  T=C02«  r  J  ) - j 

CALL  HTRT  (T,A,h(2)>  I 

HI2)=H(2>*TT/WT(2)  +96290.  *l.8/WT'2)  j 

00  1 3  J=l ,t  I 

13  AtJ)=H2R(J)  } 

CALL  HTRT  ( T , A, HI  3 ) )  i 

- Hm=wm*TT7Vff(3>  +  2077  ,"8'  *  1 .  H  /  W 1  I  3  I -  i 

DO  14  J=l,6 

14  A  I J ) =XN2R ( J  ) 

CALL  HTRT  ( T . A, hi  4  >  > 

H(4)=H(4)*TT/WT(4)  +2072.3  *l.R/WT(4) 

00  16  J=l ,6 

13  ATJ)gn?TTJ-) - 

CALL  HTRT  ( T  »  A , H ( 5  )  ) 

HI3>=n'S)*TT/WT(3>  +2074.7  *1.R/WT<5> 
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DO  16  J« l .6 

16  MJ>=H2nB<J> 

CAUL  HTRT  <T,A,H46>> 

W7S>»HI61*TTjm  (6;  +6U165.7  H.B/WHM 
HV2-H46) 

17  CONTINUE 
HNC2  -0.0 
00  18  1-1,5 

18  HNC2-HNC2  C(I>*H(D 

— ffNC-<r»HNu2;Tr.ir-c(bn - 

RETURN 

- ffNTJ — - - — - 


SUBROUTINE  SUMJT (C»H,CN»ENTH) 

twpcrcn  rtbt.*h  iA-H-,tmu - 

DIMENSION  C  4 l ) • H 1 1 ) 

SMM=0, 

no  i  1*1,5 

l  SUM-SUM  +  C I  I  )*H|  1  ) 

ENTm=SUM/cn 

RETURN - - 

END 


AEOCTR-72  89 


MAIN  PROGRAM  FOR  EXHAUST  GAS  SPRAY  COOLER 
WITH  ZERO  BLOCKAGE 

IMFUCIT  REA|.*8IA-H,0-2  > 

_ common  aa  _  _ _ 

COMMON  "  ROL  .VW.GW.GC.RV .C J,WNCA0,B<4)  ,lstA 

COMMON  |  SET.KKK.ROUNT _ _ _ _ _ 

COMMON/  '.OW/  MICH 

0|MENS10_N__£l(  10>,VL(  10>.TS(_10).0a0),YI45).  YP(45).STA< 10) 

103  FORMAT flMl » • BEG  I N  ST  A  T ION*  » 1 4  « 1 2X  » 1 2HRUN  NUMBER  ,  2A9  ) 

M»0 

200  M«M+1 

B A* . 2  5  _ _ _  _  _ 

M  1Ch*0  .  ■ 

J»EAO(5,5l015)AtR_l»A<.P2 _ _ _ _ _ 

<31015  FORMAT (  2  A9  > 

_ REAP  LS .102)  NSTA 

KEN  =  0 

_ _ _JiSl=  NSTAtX.  _ _ _ . _ _ _ 

REAO(5,10l>l$TAU>,l  =  l.NSl> 

_ READ(5»10Q)  CV.VG.IG.P  .  _ _ 

READ)  3. 1 01  >  IF'.  <1  )  .  1  =  1  .N$TA> 

_ REAP  (6.101  XVI  t I >. 1*1  .AISTA) _ 

READ! 101 H  TS( 1  I . 1*1 .NSTA ) 

RE  AO (5. 101) <D( I», 1  =  1. NSTA) 

RFAOI  S,  101  HB(I>.1*1.4> 

100  FORMAT (4E10. 2 > 

101  FORMAT!  7E10.2) 

_ LQ2 _ fcjJRM.11.1  12) _ 

CPI.  =  1  .0 
CJ=778.0 
GC  =  32.2 

1SET  =  <5  .  . . . 

Gw=  2  9.0 

_ lAkl^.ia.0 _ 

R V*  1.9R6*  77B.0 

ROL=62 .4  ...... 

TREE  =  <340.0 
ISTA=1 
NE  0=  R 

_ C _ THE  ARRAYS  1 1 S 6 0  IN  D1FFE  ARE  NOW  SET  UP _ 

Y< 1 ) sCV 
Y 1  2  >  =  V  0 
YI 3  >  =  TG 
Y I  4  )  =  P 
Y<  S  )  =E(_  (  1  ) 

Y  (  (S  )  =VL  (  1  )  _ _  . 

V<  7  >  =  T  S  (  l  ) 

_ yiB).-.  oil ) - 

WNCAO  =  Y(2)  * ( 1 ,0-Y ( 1 ) >*Y ( 4 ) / < Y ( 3 >*RV* < ( 1 -Y ( 1 ) ) /GH+Y ( l ( / VW )  > 

X  =  0.0 

WR  I  TF  ((S.  103  )  ISTA,  ALP1  ,  A'.P2 
3  CONTINIIF 

nx= .oooi 

_ KQlINT  =0 _ 

00  1  1=  1.499999 

C  A 1. 1.  0|  FFF  (  X,Y,YP,f)X,KR<,NFO,I.KFN) 

|P ( x  .OE  .STA (  1STA.1 ) )G0  TO  2 
l  CONTINUE 

,IJJ=  JJJM 
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2 _ CONT  I  NUP _ 

C  WE  WILL  NOW  PLOT  THE  STATION  JUST  FtMSHtO 

I  ST  A=  I  S  T  A  ♦  1 _  _ 

I S6  T  =9 

Iff ISTA  .GT.  NST A )GO  T0_5 
WR|TEf6.103)ISTA',ALPl,ALP2 

C _ SfT  IIP  AS  ft  A  Y  S  FOR  01  PEE _ 

VINEO+1 >*PL ( I STA ) 

V ( Nf  0  +  2 >  =  VL ( I  ST  A ) 

Y ( NFQ+3  >  =  t  S  ( ISTA> 

Y<NfQ+4>=ni ISTA> 

NEO  =NE0+4 

_ KOUNT  =  0 _ 

r.n  Tn  3 

5  IFIM-1 >200.222.222  _ 

222  STOP 

F  NO 


SO  AKOl.'T  INE  r>IFFE(X,Y,YP,r>X,<K,NFO,  f.KEN) 
IMPLICIT  PFAL'-Fl  A-h.O-Z  > 

01  MANSION  Y(50>.YP(50>.7<50>.ZP(50),ZNf50> 
CALL  YFMNC ( X, Y, YP, 1 ,KEN > 

00  1  !  =  l  .NEO 

1  2 (  I  ) = Y (  I  >  +0X  * YP I  I  > 

X=X+OX 

nx2= . sn+o*nx 
no  s 

CALL  YFIINC  (  X  ,1  .2P..I.K6N) 

K=0 

no  AO  j  =  1  . NF  0 

ZN  (  I  >  =  Y  (  I  ) +0X2 ‘  ( VP  I  I  >  +  Z  P (  I  I> 

]F(0AAS(ZN<I)-Z<1>  > -l .0-0 5*0 A AS (ZN(I>>)4,4,3 

3  K  =  l 


4  K«  =  J 

Z (I  >  =ZN (  I  > 

40  CONTINUE 
1 F ( * >5, *.5 

5  CONTINUE 

WR|TP(A,99XZN(|) ,1  =  1  ,nSO  > 

WR  I  TF  (  is  ,  «0  >  K  K 


- ^  fnPMAT(T5> - 

99  F  OR  w  A  T (4F20.10)  _ 

STOP 

6  00  7  1=1, NEO 

7  Y I  I  ) =Z  <  1  > 

_ I  F  (  J  .  C.F  .  3  .AMO,  j  ,LF.  5  )  GO  TO  1212 

IF(J  . L  T .  3  >00  TO  2020 
ox= . s«nx 
on  To  i?i2 
2020  O x  =  2 . 0 *0X 

I  c  (  Ox  .GT..01 >0X=.01 

1  2  1  2_C  ON  T  |  Ntlf _ 

KFTOPN 
L  NO 
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SUBROUTINE  YFljNC  <  X  ,Y  ,YP,K  ,KEN  > 

_ _ I_MPL  I  CJ_T_RE_AL  *6  (  A  *H, O-Z  ) 

COMMON  BA 

_ COMMON _ ROL . VW.GW .PC  ,R V.C J ,  HF'C  AO,  B(  4  )  ,  I  S  T  A _ 

COMMON  1 $6  T ,KKK ,  KOUNT 

_ 0 1  MENS  I  ON  V I  50 ) <  VP I  50 > , A5 (10 > , A6 < 10 ) , A7 I  10 > . ABI 10 ) , A • 0 <  » 0 > , X VS < 1 0 

l ) ,XHL  < 10) ,XHFGTi6>Y8ARH< 10) ,0M< 10  >,  XKI  10) , CD!  10) 

_  2 . XMV ( 10) «R6)  10  > ,SCF ( IQ) .ABNI 10  > ,PVF (10). FNU ( 10 ) 

pie  =3.i4i6 

_ SUM1=Q, _ 

S()M2  =  0. 

SUM3=0. _ _  _ 

SUM4=0. 

_  SUM5=0,  _  _ _ .  . 

CALL  ROGeeiROG,Y(l).Y(2).X) 

_ 00  1  (  =  1  .  I s ta _ 

J=  B*4*(  1-1 1 

_ A5(l>  =  PI E*ROL*Y( J )**2  /2.0  _  _ 

CALL  SUB 1 (  y'(  1  >  ,Y<  3  )  ,Y(4>  ,V(  J-l  >,XVS(  I  ) . XV.XK (  I  ) , Y IJ ) , DM (  I) , CO  I  1  )  . R 

_ l ARH ( |  >  ,  Y(  J-2  I.  Y(2_>.RE ( 1  ) f  SCF( I > , ABN ( !  > , PvF ( ! ) . FNIM  !  )  ) 

YP(J)=  XKI  I  ) *  PIE  *Y(J)**2  *1 XV-XVS(l)  )»VW/( A5(  I ) *  Y I J-2 )*( 1 -XVSI  I  ) 

_ LL2 _ 

AM  |  !=  YPI  J  ) 

YPI  J-2  )  =  (  I  PI  e/8.0  >*ROG*<  V  I  J  )  *  *2  )*f.OI  I  ) s  D  AR  S I  Y  I  2  )  — Y  I  J-2  )  )  =  (  Yl  2  )  — 

1 Y I J-2 ) ) - 1  Y(J-2)-Y(2))*Y(j-2)*AS(  ]  >*A6I  I ) ) /  I  DM  (  t  j  *  V I J-2  )  ) 

A  7  I  |  )  =  YP( J-2) 

YP( J-3 )=  AS  I  I )  *  A6  I  I  )  *  Y  ( J-3 )/DM(  I  > 

_ A  B ( [  )  =  YP( J-3) _ 

SUM l  =  SUM  1* AS (  |  ) *A6 I  I  >*Y ( J-3  > /OMI  I  ) 

call  hl  I  i  y  (  j-  i  )  •  xhl  ii).  xmfg  1 1  >,  xmv  im 

YP(J-l)  =  BARHII)*  PIE  *  Y(J)*»2  *(  Y(  31-YI  J-l  )  )/  (0V(I)=»Y(J-2))- 
CXML  (I  )*AS(  I  )«A6  I  I  )/OMI  |  )  +  A5(  I  )*A6(  I  )*XHV(  I  ) /l)M(  I) 

A’, 01  I  )=  YPI  J-l  ) 

_ SUM?  =  SUM2  ♦  Yl J-2 )»A8(  I  )♦  YIJ  — 3>»A7II) _ 

SUM3  =  Sum 3  *ah ( I ) *( XHL I  I )  +  Y ( j-2  >**2  /I 2.0«GC*CJ  )  ) 

S I IM4  =  SUMS  ♦  Y ( J- 3  )  a  Y I J-2 ) «A7(  I  )/ (GC*C J ) 

SUMS  =  SUMS  +  Y I J-3)*A10( 1 ) 

1  CONTINUE 

YPI l >  =  —  Sum 1  * | I-Yl t : )  =  *2 

_ A9=YP|  t  ) _ 

CALL  AAA3IX.YI1  )  «  Y  (  2  )  ,  A  3  ,  A  A  0  ,  D  A  A  0  <  R  OG  ) 

CALL  OERT (  ROG.  Yl  1  ) . Y (2  > , Y ( 3) .Y I  4 ) , A9, SUM?, SUM3, SUM4, SU“S, AO, A3. A  1 

Cl. A12, A13, A14, R, CP V, CPA. A?. A4, A AO. A  * 

VP  (3)  =(  AH*A14-A12*A|  3  )  /  I  A1  l*ICPV*Y(  I  )/(  l-YI  l  )  )  .CP  A ) -A  1 ?* AO*ROG* R 
C  ) 

A  1 S=  YPI  3 ) 

YPI?)  =  I  A  1  3-A0«R0G*K*A IS ) /A  1  1 

_ a i 6=yp  1 2  ) _ ; _ 

CALL  0ERPIR0G,V(1),Y(?),YI3).R.A2,A3,A4,A9,A1S,A16,DP,AI 

YP I  A )  =  OP 

IFIK.GT.  1  )  GO  TO  3 

1FIX  ,EL>.  O.OIGO  TO  9999 

IF  (  X  .L  T .  BA )G0  TO  3 

_ fiAsfiA*  .25 _ 

9999  WR  I  T  F I  6 ,4  >  X 

WR I TE I M  R 7654  )KEN 
B  7654  FORMAT  I  I  30 ' 

MM-  uli/u.o-mn 
A  FORMAT  I  1H0 , • x  =  •  , F 1 6 . 1 0 ) 

_ 9R1TE(l,S>  AAO.WH 
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5  FORMAT  (  1HO  t  *  A/AO®  *  ,G16.6,5X,  • W V/WNC* • . G16 . 6 ) 

_ WR I  TE  (  6 , 6  )  Y(1  >.Y(2).YI3)  ,Y(A)  _  _____  __ 

6  FORMAT!  1HO,  'CV»* ,G20.9,5X , •  VG= • ,G20. 9 , 5 X , ' TG» ' , G20 . 9 , 5X , • P« • ,G20. 9 

00  10  !»  I.ISTA  " 

WRITE  (6.R)  YU-3!  ,Y(J-2>,Y(J-1),YIJ> 

8  _  FORMAT  (IHO.'F  L  « •  ,  G20.  9 , 5X, >VL  *  %G20 . 9 ,  5X,,»  T$» ' ,020. 9, 5X ,  'O*1 .G20.9 

l>  . . 

_ WR|TE(6jl9>  YPl J-3> . YP f J-2 )  ,  YP< J-l  ), YP( J  ) _  ' _ ' _ _ 

WRITE! 6 ,9R9 )CD I ! I « RARH ( I ) »RE  ! I  ) ,  SCF !  !  ) . 

QB9  FORMAT!  lHQ,  <CD=,,016.6.5X,‘BARH='.G16.6.5X,’RFgt,G16.6.3X,»SCF»l, 
1016.6) 

_ WRITE (6,R99>ARN! ! >,PVFI I ) ,FN!M  I  )  _ _ 

R99  FORMAT! IMO , 1 AAN=  • ,G16.6 ,5X , •  P  VP~*  *  »G16 . 6  «  5X  , •  FNti* ' ,016.6) 

10  CONTINUE 

9  ‘  FORMAT!  1H0. 'FLP**  ,G16. 6. 5X  ,»VLP=*  ,G'16.6,5'X,"*  TSP-"'VgTo.'6‘.  5X, 'OPs*  ,G 

_ 1  16. A) _ _ 

WR I T  F ! 6 , 7  )  YP!  1  ) , YP( 2 )  ,YP ( 3) . YP( A ) 

7  FORMAT ( 1H0, *CVP=* ,G'  6,5x, *VGP= • ,016.6, 5X, * TGPe • , G 16 . 6, 5X , *PP*' ,0 

116.6) 

1001  CONTINUE 
1000  CONTINUE 
3  CONTINUE 

1 - FTiRRxrriH-;7d4-.i  > - ; - 

KEN=KEN+1 

RETURN 

FNO 


5iiRR0wT  J NE  SuRl (CV.TG.P.TS.X VS. XV, XK,0, DM. CD. RH.VL.VG.RF.SC F.ARN. 

1 P VF  «  FNU ) 

TMP1. 1 C  IT"  RF  AL  *R  I  A-H  ,n-7  I - 

C  OMMON  AA 

COMMON  ROL.VW.OW.OC.RV.C J.WNCAO.RIA)  ,l$TA 

COMMON  ISFT.KRK 
XV=lf.V/VW)/!CV/VW  ♦(  l-CV  >/GW  ) 

x?=-9.06* ( ! -.S69fe*TS*  .0R39E-A  »TS*»2  ♦.0927E-7  * T S** 3* 1 352 . 3 J / T S- 
C I  ,AS?6744 I  ■ 

X  I =  67?. 0/ T  9 


PVTS  =  ?UT.O  •Xl*«5.19«f)EXP(  xT> 
xvl  =  PVTS/P 

X  V  E  =(  XVS +  XV  l/?.0 

XMi  XV  F.‘  VW  ,(  I-XVE  I  »GW 

TE!  =(  TS->  Tf, )/?  .0 

FMV  s(10.6E-7  *OSORT ( TF |  )  ) / ( 1*1 53R.O/TF 1  ) 

- - 

FMNC=  7 . 5E  —  7*  O'iOR  T  (  T  E  I  )/ll*  216.0/TFI) 

EM  =  X  VE  *FMV  ♦ ( 1 -X VF  ) »FMNC 

OAR  =  1 .5375/P)«(TE1/491.0)**2.33A 

RGF  =  P*  XM / ( R  v» TF  |  ) 

SCE-EM/  (  RI)E  #|)AB  ! 

RE  =0A5S(  V(.-VL  )»ROE»0/FM 

ARN*  ? . 0+0. 6* n SORT (RE )  *  SCE»*  0.3333 
*K  =  AHN*FM/ | 0*SCE« XM > 

OM  =  n«-»x  *  3 . 1 A  1 6  "ROL/6.0 

C.o  =  ?..0/Rc«(  1.0*.15*RF-x*0.6R7) 

1E(  T  F  |  .0T.17C10.0  .AND.  T  E I  .LT.A500.0J  00  TO  221 
I  E  (  T  E  T  .It.  RQQ.O  .OR .  TE  »  . G  T . a  500  )  GO  TO  222 
CRVI-  =.A3^A  +.167ftE-A*TE|  *0. 2  7R  1E-7«TF  1**2. 0 

CPNC?  =  0 . ?  3 1 R  +  O.IOAOE-a  #TFI  ♦  0.7166F-R  *TE|*«2 
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an  m  102 

221  CPVF  =  _,3_3J9_*0,  14386-3  “TFI  -0.13126-7  *TFI**2.0 
CPNCE=  0.2;'i4  ♦.  0.3521£-4»TF  1-  0.377AE-8  *TFI**2 

_ an  to  102 _ 

222  WRlTFIft.lOll 

101  EORMATUH  ,»,TEMP-F|-|S  OUT  OF  ft  A  NG  EJ  )  _  ■_ 

WRITE (6, 122)  TG.TS 

122  .  FORMAT ( 1  HQ . 1 TG“ 1 .620.6.  5X , « TS“ » . 620. 6 ) _ 

STOP 

1 02  CONTINUE _ 

FKV  =  0 . 4  32  “E  MV 

FKMF  -  Q.2S7  *10.115+  5.17»CPNCF>»  FMKlC  _ 

E<  =  XVF»FKV  + ( 1— XVF ) *  FKNC 

CPF  =  XV  F-MVW/XM  )«CPVF  » ( 1-XVF ) »CPMCF*OW/XH  _ 
PVF  =  CPF  “EM/FK ' 

FNU  =  2* . 00  «OSQRT(RF)  ♦  PVF»“  0.3333 _ 

BH  =FNU  *F  M  *CPF/(0*PVF) 

RETURN 

FNO 


SUBROUTINE  ROGFE (ROG.CV.VG.X  ) 

IMPLICIT  «6AL*8( A-H.n-Z ) 

COMMON  BA _ _ _ _ _ 

COMMON  ROL.VW.aw.GC.RV.CJ.WNCAO.BIA)  .ISTA 

AAO  =B(1)  +8<2>»X  ♦  B  ( 3 ) “X  *  *2 • 0+  B(4)*X**3.0 
ROG  =WNCAO/(  (  l-CV  )*\/G*AAO> 

RETURN 

ENO 


SUBROUTINE  HLIITS.XH1  .XhEG.XHV) _ 

IMPLICIT  «EAL*8( A-H.0-2 ) 

XHL  =  T  S  -540.0 

XHFf,  =  - . 5896  *  T  S  +  .08396-4  *TS**2  ♦  0.0927E-7  * T S*“ 3 . 0 ♦ l 35 2 . 3 
xhv=xhL*XhFO 

RE  TURN 
ENO 


SUBROUTINE  AAA3(X,CV,Va.A3.AAO.OAAO,ROG) 

IMPLICIT  RFAL+BI A-h.0-7 > 

COMMON  BA 

COMMON  ROl,  1  WW  tGW  ,G0  ,R  V  ,C  J  j  WNC  AO,  B(  4  )  ,ISTA 

AAO  =  Bill  ♦  B(2)“X  +R(3>*X**2.0  *R I  4 )« X““3.0 
OA AO  =  R<?)  *2.0“  B (  3  t  *  X  *3.0“B(4I“X**2.0_ 

A3  =  VG  *I1-CV)»  OA  AO/  1 VG»I  1 -C V ) » A AO > *“ 2 

SETUKN - 

ENO 


Subroutine  ofrt<  rog.c v , vg.tg, p, A9, sum2 . SUM3, siim4, sums, ao. a 3, a 1 1 , 
1A12,A13.A14,R,CPV,CPA.A2.A4,AA0,A) 

IMPLICIT  RE AL * M A-H , 0-7  I 
COMMON  BA 

COMMON _ RPL,VW,0WiGC>BV.CJ.WNCA0.B(4)  ,ISTA _ 

COMMON  |  Sc  T  , <<R 
IMVG  .  L  T  .  0.0  )  G  0  TO  70 
Rt  P/(ROf,*  TG  ) 
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AO  _=  GC»(  l  -C  V  >  /ROG  _ _  _  _ 

A*  *  WNCA0/(VG**2*!1-CV>  *AAO> 

All  s  ( 1.0-CV>»VO-TG»R»A0»A* _ 

A12=  C1*CV/<1-CV> >*VG/(GC*CJ  > 

Al  =  GW  *C V  ♦VW»(1-CV>  ____ _  _  _ _ 

A  =  WNC  AO/  I VG* ( 1 -C V  j  **2~*AA0  > 

XV*  GW  *C  V/ A  1 _ _ _  _  _  _ 

A2  =  RV»VW*GW*!GW-VW  >/ I  Al**2'*<  XV»VW*<  i-XV)*GW)**2  ) 

A  1 3  =  A0»WNCA0»TG*R»A3-A9»( VG»»2 +ROG* TG» A?«AO+TG*R*A  »A0)  -VG«(1-C 
CV  >**2*$()M2 

CALL  hFFV  < tg.xhv.cpv.cpa  )  __ 

414  ®-A9*'(  XHV*VG*»2ro/ l2*GC*C.m/ l  i-CV)**  2  - SUM 3-  SUM*-'$(jM5 

R  6  TURN 

20  WR|TF!6,3l  )VG 

3i  format<f?o.io> 

STOP 

FNO 


SUHROU  T I NF  hFFVI TG.XHV.CPV.CPA ( 

IMPLICIT  REAL*B<A-H,0-Z > 

COMMON  8A  "  . 

COMMON  ROL. VW.GW.GC ,Rv,C J, WNCAO.BI 4)  ,!STA 

COMMON  I  SET  ,KKK 

IF(  TG  .OF.  1  700.0  .ANO.  TO  .LF. <-500>Gn  TO  222 
IF<  TG  .i,T.*00.0  .OR.  TO  .GT.4S00)  GO  TO  221 . 

X H V  =  ,*30*“TG  ♦. OR39E-**TG**?. 0+0. 0927F-7«TG**3. 0-236. 3161+10*2.9 
CPV  =  .430*  ♦  . 1678F-**TG+.27R1F-7*TG**2.0 
CPA  =  .7318+  .10*0F-**TG+.7166F-R  »TG**2.0 
rn  Tn  10Q 

22?  XHV  =  .3319  *  TG  ♦  .0719E-3°TG**2.0-«0*373E-7*TG**3-1R5. 381+10*7,9 
CPV  =  .3319  ♦  . 1 *3hE-3*  TG-  0. 1 3 1 2E-7* TG** 2 . 0 
CPA  =  .221*  +0. 352 IF-**TG  -0. 37 76F-8* TG** 2. 0 
GO  TO  IO0 
221  WR 1 TF | 6. 99  ) 

7^5  pormat ( ih  ,*7EmA-g  Is  hot  op  range1 > 

WR I TF ( 6, 1  )  TG 

l  FORMAT! IMO.'TG*' ,F16. 6) 

STOP 

ion  cnNTiNne 

RETURN 

FNO 


SU8 ROUT  INF  DERPIROG.CV.VG, TG,R,A2,A3.A*,A9,A15,A16,DP,A) 

TBPr J C  T  T  REAL *A ( A-H  ,0-1  ) 

COMMON  RA 

COMMON  ROL  .VW.Gw.GC.RV.C 'j.WNCAO,  8!  *")  .  IST  A 

()P=  ROG*R  +  A l 5  +(R0G*TG*A2  +  TG*R*A  I  * A9- TG* R *  A** A 1 6-WNC A G > A 3*TG*R 

RF  TURN  "  . . 

F  NO 


